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i. introduction and summary 

As part of its V/STOL research program, NASA intends to conduct flight 
investigations of the stability, control and handling qualities of highly 
augmented V/STOL aircraft. Specific plans include the flight tests of a 
YAV-8B aircraft. Figures 1-1 and 1-2, modified to include an advanced 
avionics and flight control system for improved flying qualities and perform- 
ance. 


As an initial phase to this program, NASA will conduct flight tests of 
the YAV-8B vehicle in order to extract aerodynamic and propulsion character- 
istics, update existing simulation models, validate handling qualities and 
design criteria, and to improve V/STOL flight test techniques. This program 
will also include tests using a static test stand located at the Dryden 
Flight Research Center, where flight tests of the YAV-8B will take place. 

In order to perform high quality parameter estimation and analysis of 
the YAV-8B characteristics, it is necessary to construct mathematical models 
of varying complexity and linearity from existing wind tunnel and flight 
test data. 

The McDonnell Aircraft Company (MCAIR) recently completed a V/STOL simu- 
lation and modeling study under contract to NASA Dryden. This study defined 
and documented detailed mathematical models of varying complexity representa- 
tive of the YAV-8B aircraft. These models will be used by NASA in parameter 
estimation and in linear analysis computer programs while investigating 
YAV-8B aircraft handling qualities. Both a six degree of freedom nonlinear 
model and a linearized three degree of freedom longitudinal and lateral 
directional model were developed. 

The nonlinear model is based on the mathematical model used on the MCAIR 
YAV-8B manned flight simulator. This simulator model has undergone periodic 
updating based on the results of approximately 360 YAV-8B flights and 8000 
hours of wind tunnel testing. Qualified YAV-8B flight test pilots have 
commented that the handling qualities characteristics of the simulator are 
quite representative of the real aircraft. These comments are validated 
herein by comparing data from both static and dynamic flight test maneuvers 
to the same obtained using the nonlinear program. 

The linearized mathematical model uses stability derivatives and is 
formatted exactly as the models traditionally used in conventional flight 
dynamic analysis. Aircraft characteristics were predicted using this 
linearized model and compared to both the flight data and the nonlinear 
predictions. To document the aircraft characteristics throughout the flight 
envelope trim conditions and stability derivatives are provided for 24 flight 
conditions. 

A FORTRAN batch simulation of the nonlinear model nas been produced. 
Documentation for this simulation consists of a description of the software 
including top level flow charts, program structure, subroutine interfaces, 
modeling equations, data format, a user's guide, source listings and plots 
of the over 17,000 aerodynamic and propulsion data points used in the model. 
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This report is divided into two volumes. Volume I contains the descrip- 
tion of the aircraft, documentation of the nonlinear and linear mathematical 
models, stability derivatives, the comparisons of predicted and actual flight 
test data which validate the nonlinear model and a discussion of models 
appropriate for use in parameter estimation programs. Volume II contains the 
source listings for the nonlinear program and plots of the aerodynamic and 
propulsion data used in the nonlinear program. 
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2. 1 GENERAL ARRANGEMENT 

The YAV-8B is a single seat transonic light attack V/STOL aircraft 
powered by a single YF402-RR-404 turbo-fan engine. The YAV-8B aircraft is 
an advanced version of the AV-8A aircraft with greatly expanded and improved 
capabilities. It retains the characteristic appearance of the AV-8A while 
incorporating an improved inlet design, a larger wing with an advanced tech- 
nology airfoil, zero degree scarf of the forward nozzles, and under-fuselage 
modifications called Lift Improvement Devices (.LIDs) . A three-view general 
arrangement drawing of the YAV-8B is shown in Figure 2.1-1. 

Conventional aerodynamic controls are utilized in wingbome flight and 
engine bleed air reaction controls are used in jetborne flight, with both 
systems operative during transition modes. 

An all-movable tailplcne and front and rear pitch reaction control 
system (RCS) valves, blowing downward, provide longitudinal control. Lateral 
control is provided by outboard ailerons and ECS valves, blowing down or up 
and down for large lateral inputs. Directional control is provided by a 
conventional rudder and a yaw RCS valve in the aft fuselage tail cone, 
blowing sideways, left or right. 

In the high lift mode for VTO and STO, the ailerons are drooped 15 
degrees and the flaps are interconnected with the engine nozzle control. 

This - ’ 1 ;s maximum longitudinal acceleration and good longitudinal control 
mar o be obtained during ground runs and transition to wlngborne flight. 

Shot adings and vertical landings can be made with either the mid 25° or 

the 61.7 flap setting. Conventional landings are made with the 25° 

flap position. 

A maneuvering flap, which can be positioned at any deflection from zero 
to 25 , i <s provided for increased maneuverability throughout the flight 
envelope. 

The LIDs consist of a retractable fuselage fence located at the forward 
end of the gun pods and two fixed strakes on the gun pods. Larger strakes 
are provided for attachment to the fuselage when gun pods are not carried. 

The YAV-8B has provisions for external store loadings on six wing sta- 
tions and the fuselage centerline. In addition, guns in two pod can be 
attached to the fuselage. Recent flight testing of the YAV-8B has been 
performed using a removable Leading Edge. Root Extension (LERX) . This device, 
designed to increase the wlngborne maneuverability, has only a small effect 
on the V/STOL flight characteristics. Since considerably more flight data 
exists on the aircraft without the LERX this report assumes the aircraft to 
be configured without the LERX in the Basic configuration of (.5) pylons and 
( 2 ) gun pods. 
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2.2 DIMENSIONAL DATA 

Details cf YAV-8B physical characteristics are presented in this section. 


2.2.0 WEIGHT AND FUEL CAPACITY 

Operating Weight, Flight Test Aircraft - Lb 13847 

Internal Fuel - Lb 7185 

2.2.1 WING DIMENSIONAL DATA 

Area (theorectical) - Ft* 230.0 

Span (projected) - Ft 30.33 

Aspect Ratio 4 . 0 

Taper Ratio 0.3 

Chords (projected) : 

Root (theoretical) - In 139.99 

Tip (theoretical) - In 42.0 

Mean Aerodynamic - In 99.79 

Thickness Ratio (t/'c) (FIGURE 2.2-1) 

Root - % 11.5 

Tip - % 7.6 

Incidence (w3 th respect to FRL) : 

Root (exposed theoretical) - Deg. 3.0 

Tip - Deg. -5.0 

Sweepback (projected): 

Leading Edge - Deg. 36.0 

1/4 Chord Line - Deg. 30.62 

Dihedral - Deg. -11 

T/fist (FIGURE 2.2-2) - Deg. -8 

Airfoil (FIGURES 2.2-3 and 2.2-4) Modified 


Supercritical 

2.2.2 HORIZONTAL TAIL DIMENSIONAL DATA 


Area (projected) - Ft^ 47.54 

Span (projected) - Ft 13.92 

Aspect Ratio 4.079 

Taper Ratio 0.201 
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Sweepbac k ( pro j ec ted ) : 


Leading Edge (theoretical) * Deg. 39.81 
1/4 Chord Line (theoretical) - Deg. 33.91 
Trailing Edge - Deg. 10.77 

Dihedral - Deg. -15.84 


Chords: 

Root (centerline of airplane) - In. 
Tip (theoretical) - In. 

Mean Aerodynamic chord - In. 

Pivot Line Location: % MAC-- 


(F.S . 559.05 and W.L. 22.35) 22.42 

Airfoil Section: RSA Symmetric 

Thickness Ratio (t/c): 

Root (centerline of airplane) .07 

Tip (theoretical) .07 

Tail Length (from 25% wing MAC to 25% tail MAC) - In. 199.01 

Horizontal Tail Volume Ratio .412 


Deflection (with respect to water line plane) 
(Includes +1.5° Autostab) 

Maximum Mb se down - Deg . 

Maximum Nose up - Deg. 

.2.3 VERTICAL TAIL DIMENSIONAL DATA 

2 

Area (theoretical above W.L. 126.50) - Ft 
Span (W.L. 126.50 to W.L. 194.00) - Ft 
Aspect Ratio 
Taper Ratio 

Chord : 


Root (W.L. 126.50) - In. 86.93 

Tip (theoretical) (W.L. 194.00) - In 23.3 

Mean Aerodynamic Chord (W.L. 153.76) - In 61.24 

Airfoil Section: HSA Symmetric 

Thickness Ratio (t/c) : 

Root (W.L. 126.50) .082 

Tip (theoretical) (W.L. 194.00) .052 


25.83 

5.63 

1.23 

.268 


-11.75 

+12.75 


67.39 

13.56 

46.44 
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Sweep back: 


original pag:: u 

° F P °0R Q'JAlITY 


Leading Edge - Deg. 

1/4 Chord - Deg. 

Trailing Edge - Deg. 

Tall Length (from 25X wing MAC to 25X tail 
MAC) - In 

Vertical Tail Volume Ratio 
2.2.4 FUSELAGE DIMENSIONAL DATA 


Length (maximum) (parallel to water line) - Ft 
Width (maximum) - Ft 

Depth (maximum) - Ft 2 

Maximum Frontal Area (fuselage alone) - Ft ^ 
Maximum Frontal Area (Including inlets) - Ft 
Thrust Axis (zero nozzle deflection) (with 
respect to W.L.) - Deg. 

2.2.5 FLAP DATA 

Area (per side) 

Span (per side) 

Deflection: 

STO Itode 

a T > 50° 

K .< 25° __ 0 . 

NORMAL Mode - Leg. 

2.2.6 AILERON DATA 

2 

Area (per side) - Ft 
Span (per side) - In 
Deflection (including + 2° Autostab) 

Undrooped - Deg. 

With 15° droop - Deg. 

NOTE (Lateral stick travel restricted to 751 
of full deflection above 250 KIAS) 

2.2.7 RUDDER DATA 

2 

Area - Ft 
Span - In 

Hinge Line Location - Z W.L. Chord 

Soot (W.L. 126.50) 

Tip (W.L. 187.25) 

Sweepback of Hinge Line - Deg 

Deflection - Deg 


- Ft 

- In 


- Deg. 

- Dee . 


47.36 

40.37 
8.15 


187.99 

.058 


42.89 

8.00 

5.60 

19.2 

31.5 

1.50 


15.49 

64.54 


61.7 

25 

0-25 


6.19 

58.9 

+12/ -27 
+27.2/-10.4 


5.27 

60.75 


11.87 

33.68 

12.73 
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.2.8 REACTION CONTROL SYSTEM 

Front Fuselage Pitch Nozzle (thrusting upward) 


Position - In F.S. 137.48 

W.L. 70.36 

B.L. 0.00 

Thrust Incidence (with respect t j f.L.) - Deg 97.63 

(with respect t) B.L.) - Deg 0.00 


Rear Fuselage Pitch Nozzle (thrust in} upward) 


Position - In F.S. 

W.L. 

B.L. 

Thrust Incidence (with respect to l.L.) - Deg 
(with respect to 3.L.) - Deg 

Rear Fuselage Yaw Nozzles (thrusting port and starboard) 

Position - In F.S. 

W.L. 

B.L. 

Thrust Incidence (with respect to B.L.) - Deg 
(with respect to W.L.) - Deg 

Wing (port and starboard) roll nozzles 

(thrusting upward) 


605.51 

114.63 

0.00 

82.00 

0.00 


599.00 

120.50 

0.00 

+90.00 

0.00 


Position - In F.S. 432.45 

W.L. 83.89 

B.L. 177.84 

Thrust Incidence (with respect to V.T..) - Deg 83.00 

(with respect to B.L.) - Deg Outbd 5.00 


Wing (port and starboard) roll nozzles 

(thrusting downward) 


Position - In F.S. 432.45 

W.L. 83.89 

B.L. 177.84 

Thrust Incidence (with respect to W.L.) - 263.00 

(with respect to 3.L.) Deg Inbd 5.00 
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Reaction Control System moments versus control surface deflection or 
stick, position 

0 % of -nav-innim moment for stab ila tor @ 2° 

100% of maximum moment for stabilator 3 + 10° and -3° 

0% of maximum moments for neutral lateral stick 
100% of maximum downblowing moments fo: 55% (2.31 in) of maximum 
lateral stick deflection 

0% of maximum upblowing moments for neutral to 55% (2.31 in) of 
maximum lateral stick deflection 
100% of maximum upblowing moments for 80% (3.365 in) lateral 
stick deflection 

0% maximum moments for rudder @ 0® 

100% of max imum moments for rudder 5 +10° and -10° 



FIGURE 2.2-1 

YAV-8B SPANWISE THICKNESS DISTRIBUTION COMPARISON 
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FIGURE 2.2-2 

YAV-8B GEOMETRIC TWST DISTRIBUTION 
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FIGURE 2.2*3 

YAV-8B DESIGN STREAMWISE CAMBER DISTRIBUTION 
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FIGURE 2.2-4 

YAV-8B LEADING EDGE RADIUS DISTRIBUTION 


2.3 WEIGHT AND BALANCE DATA - Included within this section are the weight, 
inertia, and center of gravity data for the YAV-8B flight test configuration 
oi. clean aircraft with (2) gun pods, (2) inboard pylons, (2) intermediate 
pylons, and centerline rack. The moment of inertia and C.G. position data 
for this configuration are shown in Figure 2.3-1. In addition to the effect 
of pylons and gun pods, the zero fuel weight includes the pilot, oil, trapped 
fuel, and oxygen. A center of gravity envelope is presented in Figure 2.3-2. 
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1 . Takeoff gross weight 

2. Water expended 

3. Wing fuel partially expended 

4. Combat weight (40% fuel expended} 

5. Wing fuel expended 

6. Fuselage fuel partially expended 

7. Fuselage fuel expended 

8. Feed tank fuel expended 



2 4 6 8 10 12 14 16 18 


C.G. • percent c 


OMi-o ttt-n 


FIGURE 2.3-2 

YAV-8B FLIGHT TEST AIRCRAFT 
CENTER OF GRAVITY ENVELOPE 
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2.4 FLIGHT CONTROL SYSTEM 


The YAV-8B flight control system consists of conventional ailerons, 
rudder, and stabilator, with a reaction control system (RCS) acting about 
all three axes during hover and transition. The s^bilator and ailerons 
are power operated while the rudder is connected directly to the rudder 
pedals. In view of the ineffectiveness of aerodynamic controls at low 
speeds, high pressure compressor bleed air is fed to the aircraft extremi- 
ties and ejected through variable shutter valves connected to the conven- 
tional controls to provide reaction control. 

A single channel limited authority stability augmentation system (SAS) 

Is provided to facilitate control in hover and transition. It may be 
selected "on" whenever the aircraft is flying below 250 KIAS with either 
the gear or flaps extended. 

2.4.1 LONGITUDINAL CONTROL SYSTEM - Longitudinal control of the aircraft is 
maintained by a combination of a stabilator and reaction control system as 
shown in Figure 2.4-1. Control stick movement is transferred to the stabi- 
lator as shown in Figure 2.4-2. Longitudinal feel is provided by a hydraulic 
q-feel unit, augmented by nonlinear backup springs and a bobweight. The 
bobweight contributes 1.5 lb/g at the stick, generated by vertical g forces. 
The stick force gradient and dynamic pressure relationship is shown in 
Figure 2.4-3. This q-feel unit is operative above 250 KIAS and provides a 
stick stiffness increasing in proportion to dynamic pressure. Longitudinal 
trim is obtained by electrically biasing the q-feel unit which ensures full 
control authority. Trim capability Is from -4 to +7.5 degrees of stabilator 
and operates at 2.2 degrees/sec. The effect of trim range on spring force 

is shown in Figure 2.4-4. 

Low speed pitch control is by reaction control valves at the nose and 
at the rear fuselage extension, the latter being fed through the same ducting 
as the yaw valves. The front pitch shutter is linked to the control column, 
and the rear shutter is linked to the stabilator. The relation of stabilator/ 
stick position to pitch valve opening is shown in Figure 2.4-5. The pitch 
RCV valves thrust upward. 

The pitch autostabilizer has an authority of +1.5 degrees of stabilator 
and works through the forward and aft reaction control valves. The forward 
RCS SAS can be disaxigaged separately. 

2.4.2 LATERAL CONTROL SYSTEM - Conventional axlerons in combination with 
reaction controls provide lateral control as shown in Figure 2.4-6. The 
ailerons are power -operated with artificial feel derived from a nonlinear 
spring, which can be biased by an electrical trim motor. Trim 
capability is equivalent to 1/3 of lateral stick deflection and operates at 
1 deg/sec of aileron. The aileron deflection and spring force due to stick 
movement are shown In Figure 2.4-7. Lateral stick travel is restricted to 
75% of full deflection above 250 KIAS, but the pilot can override the limiter 
in an emergency. 
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FIGURE 2.4*2 

YAV-8B STABSLATOR ANGLE vs STICK POSITION 



FIGURE 2.4*3 

YAV-8B LONGITUDINAL Q - FEEL SYSTEM 

For V ca j > 250 kts 
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FIGURE 2.44 

YAV-8B EFFECT OF TRIM RANGE ON 
LONGITUDINAL SPRING FORCE 
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FIGURE 2.45 

YAV-8B LONGITUDINAL RCS VALVE OPENINGS 
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FIGURE 2.4-6 

YAV-88 LATERAL CONTROL SYSTEM 
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Roll control valves on the wing tips are linked to the ailerons and 
thrust either up or down. Lateral stick position versus roll valve opening 
is shown in Figure 2.4-8. 

A roll auto stabilizer, similar to the pitch system, is also provided 
in the stability augmentation system (SAS) and has an authority of +2° of 
aileron. 

2.4.3 DIRECTIONAL CONTROL SYSTEM - Directional control is provided by a 
combination of a conventional unpowered rudder and directional reaction 
controls as shown in Figure 2.4-5. Rudder feel is primarily aerodynamic 
but a centering spring is fitted to produce some feel at low speeds. Rudder 
deflection and spring force versus pedal deflection are shown in Figure 
2.4-10. The rudder trim has a capability of +9.5° of deflection. 



- 30 _ 20 -10 0 10 20 30 


Aileron deflection, 6 A * deg 

QP21 -0902-11 


FIGURE 2.4*8 

YAV-8B LATERAL RCS VALVE OPENINGS 
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Rudder angle deg 



0P21-0 

FIGURE 2.4-10 

YAV-8B DIRECTIONAL CONTROL SYSTEM CHARACTERISTICS 


Reaction control in yaw is by port and starboard valves in the rear 
fuselage extension, the shutter being opened by rudder movement as shown in 
Figure 2.4-11. 

The yaw autostabilizer is provided in the SAS and has a yaw reaction 
jet authority of +5 degrees of equivalent rudder. 

2.4.4 FLAPS AND DROOPED AILERONS - The YAV-3B high lift system consists of 
single slotted flaps and drooped ailerons. Operation of the flap and aileron 
droop system is controlled by three switches (flap actuation switch, 
flap mode switch and flap position switch) . The flap switches and system 
logic are presented in Figure 2.4-12 and discussed below. 

2. 4. 4.1 Flap Actuation Switch - This 3-position, toggle lever-lock, flap 
switch is located on the left aft console. Switch positions are ON (neutral), 
OFF (aft) and RESET (forward) . Selection of the ON position activates the 
entire flap system and locks the switch out of the OFF pc sition with the 
lever-lock. RESET is a momentary position which reactivates the normal 
modes of flap operation if the system has failed due to a transient electrical 
problem. Flap asymmetry must be corrected before RESETting the flap system. 
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FIGURE 2.4-11 
YAV-8B 

DIRECTIONAL RCS NOZZLE OPENINGS 



FIGURE 2.4-12 

YAV-88 FLAP SYSTEM/AILERON DROOP LOGIC 
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2. 4. 4. 2 Flap Mode Switch - This 2-position switch is located on the left 
forward console. Switch positions are DOWN and DP. The DOWN position is 
used for both takeoff and landing. The DOWN mode switch position activates 
the flap-nozzle interconnect schedule (.shown in Figure 2.4-13) if the air- 
speed is less than approximately 165 KIAS (actual airspeed limits are 175 
K1AS accelerating and 155 KIAS decelerating). The flaps move down at 45.9 
degrees per second and up at 64 degrees per second (no load rates) . The 
flap-nozzle interconnect is overridden at higher than the above-mentioned 
airspeeds, and the flaps are driven to the 25 degrees position at a 7 degree 
per second no-load rate. The maximum nozzle rotation rate is about 90 
degrees per second. 

The DOWN flap mode switch position also activates the aileron droop 
system, deflecting the ailerons down 15 degrees (at a no-load rate of 8 
degrees per second) , whenever the pilot commands the nozzles to 1 6 degrees or 
greater and the airspeed is less than approximately 165 KIAS (actual airspeed 
limits are 175 KIAS accelerating and 155 KIAS decelerating) . Aileron 
retraction is at 5.5 degrees per second no-load rate. Aileron droop is 
available at any nozzle angle while on the ground by a wei on wheels 
switch with a 3 second time delay. This system allows che.xout on the 
ground, and assures full droop at the moment of nozzle rotation during 
short takeoffs. The 3 second time delay prevents premature droop retraction 
if the aircraft becomes light on the gear during lightweight takeoffs prior 
to nozzle rotation. 



FIGURE 2.4*13 

FLAP-NOZZLE INTERCONNECT SCHEDULE 

Static Relative Angles 
Flap Mode Switch - DOWN 
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When UP, the switch energizes the flap position switch on the throttle 
handle allowing pilot control of flap position. If UP flap mode is selected 
from the DOWN mode, the flaps will retract to 0° at a 7° pe£ second no-load 
rate and the ailerons will undroop at a no-load rate of 5.5° per second. 


2.4 .4.3 Flap Position Switch - This 3-position switch on the throttle lever 
is thumb operated, moment ar.’ spring-loaded to the OFF position. Switch posi- 
tions are UP (forward) , DOWN (aft) and OFF (neutral) . With the UP flap mode 
selected, this switch allows the pilot to select any flap position from 0° 

to 25° at a 7° per second no-load rate. Flap travel continues so long as 
the pilot displaces the switch in the selected direction or until 0° or 25° 
flap travel is reached. 

2. 4. 4. 4 Flap Asymmetry - If flap position asymmetry reaches 5*, the flaps 
stop at the failed position. The flaps can be raised by using the throttle 
lever switch. No asymmetry comparison is provided in this mode. 

2.4.5 LIFT IMPROVEMENT DEVICE SYSTEM (LIDS) - The LIDS consists of a 
retractable fence located between and forward of the gun pods and external 
strakes mounted on the gun pod fairings. These strakes and fence provide 
lift improvement in ground effect during VTO operation and inhibit engine 
exhaust reingestion. Below 165 PAS , the retractable fence extends and 
retracts with the landing gear. Above 165 ITT AS, the fence automatically 
retracts. 


2. 4. 5.1 LIDS Switch - A two-position LIDS switch labeled AUTO and RETRACT is 
on the lower left side of the main instrument panel. Below 165 OAS with 
the switch in AUTO, the LIDS fence extends and retracts with the landing 
gear. Above 165 KIAS or with the switch in RETRACT, the LIDs fence retracts 
or will not extend when the landing gear is down or selected down. The 
LIDS switch should be in AUTO except for conventional takeoffs and landings. 
In emergency operation the LIDS fence extends pneumatically with the landing 
gear regardless of the LIDS switch position or airspeed. 


2.4.6 REACTION CONTROL SYSTEM (RCS) - Compressor delivery (eighth stage) 
air is available on demand for the aircraft reaction control system. The 
quantity dema n ded is a function of pilot control inputs. The master valve 
is mounted on the underside of the engine, and is geared to the gain nozzle 
rotation so that maximum airflow is provided at approximately 15 nozzle 
deflection. The bleed air is ducted to the nose (front pitch valve) , tail 
(rear pitch and yaw valves) and wing tip (roll valves) , as shown in Figure 
2.4-14. Locations of the reaction control valves are presented in Section 
2 . 2 . 8 . 
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FIGURE 2.4-14 

YAV-8B REACTION CONTROL SYSTEM 
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The Scebllicy Augmentation System (SAS) of the YAV-8B is a limited- 
authority three-axis simplex autostabilization system. The SAS is provided 
to improve controllability in hovering and transitional flight conditions . 

The SAS is controlled by the AUTO STAB MASTER and AUTO STAB ENGAGE switches 
located on the left console. A separate switch is provided for the forward 
reaction control valve (RCV) series servo. Although the axes cannot be 
selected individually, failure oi one channel does not affect operation of 
the others except that the yaw SAS is inoperative without the roll SAS. 

Position switches in the flap and landing gear linkage disengage the SAS when 
the flaps are selected UP, or the landing gear is retracted. When the aircraft 
weight is on the main gear the yaw SAS is disengaged. If the flaps are selected 

DOWN, the retraction of the landing gear does not disengage the SAS. Above 

250 KXAS a signal from the air data computer disengages the SAS. The 
authority of the system about each axis is shown in Figure 2.5-1. 



SAS AUTHORITY 

AXIS 

A CONTROL SURFACE 
DEG 

A RCS OPENING 
Z OF FULL OPEN 

_ TTf _ FORWARD RCS 

PITCH AFT RCS 

-1.5 

28 

+1.5 

19 

ROLL 

+2.0 

16 

YAW 

1 

+5.0 

50 


FIGURE 2.5-1 

STABILITY AUGMENTATION SYSTEM AUTHORITY 


Functional block diagrams for the longitudinal and lateral -directional 
stability augmentation control systems are shown in Figures 2.5-2 - 2.5-4. 
Feedback gains are summarized in Figure 2.5-5 to 2.5-7. 

Figure 2.5-2 shows the longitudinal stability augmentation control sys- 
tem. The pitch rate feedback loop is closed through the stabilizer actuator. 
The aft pitch reaction control valve is driven mechanically by the stabilizer 
actuator and exhausts downward only. A series servo actuator on the forward 
pitch RCV enhances stability augmentation whan the forward RCV is demanded 
open and the RCS is activated such as in the hover mode or during short 
takeoff and accelerating transition flight. 
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A diagram for the roll and yaw stability augmentation system is shown 
in Figure 2.5-3 and 2.5-4. It consists of a rate feedback loop for roll and 
cancelled rate damping and turn coordination for yaw. Cancelled rate damping 
consists of a yaw rate feedback processed through a high-pass filter which 
cancels the steady state signal. This feedback is combined with a lateral 
at .deration feedback in a manner to improve turn coordination at the 
transition speeds. The yaw SAS only operates through the yaw RCV; therefore, 
the yaw SAS is only effective when the RCS is activated. 

Selected gains for stability augmentation are summarized in Figures 2.5-5 
to 2.5-7. 
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FIGURE 2 6 4 

Y W*B DIRECTIONAL STABILITY AUGMENTATION SYSTEM 
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PARAMETER 

VALUE 

UNITS 

\FT GEARING 

2. 75-1. 758274*6 
+. 09038152 6 0 2 
-.01149626 « 0 3 

1.085 (6 h -2.0) 

DEG/ IN 

AFT RCV GEARING 

IN 2 /DEG 

FWD GEARING 

2.126 (6 0 - .5) 

DEG/ IN 

FWD RCS SERIES SERVO 

1/(1 + .0175S) 

- 

FWD RCV CONVERSION 

1.205 

IN 2 /DEG 



DEG/DEG/SEC 

PITCH DEADSPACE 

+. 06 

IN 

PITCH RATE STRUCTURAL FILTER (1/(1 + .0222S))* 



(1/(1 + .0448S)) 

• 

STABILAIOR ACTUATOR 

1/(1 +.025S) 

- 


FIGURE 2.5-5 


LONGITUDINAL 

STABILITY AUGMENTATION PARAMETER VALUES 
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PARAMETER 

AILERON ACTUATOR 
DOWNBLOWING GEARING 

*P 

PORT GEARING 

ROLL RATE COMPENSATION 
ROLL RATE STRUCTURAL FILTER 
STARBOARD GEARING 

UP BLOWING GEARING 


VALUES 

i/a + . oi3S) 

.778 

.346 

0 + 3.438886$ 

- . 4106176* 2 

+ .04090556* 3 

a + .781S)/a + 1.29S) 

1/a + . 0164S) 

0 - 3.438886$ 
-.4106175*2 
-. 04090556 $ 3 

.477 


FIGURE 2.5-6 

LATERAL S TABILIT Y AUGMENTATION PARAMETER VALUES 
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UNITS 

IN 2 /DEG 

DEG/DEG/SEC 

DEG/IN 


DEG/IN 

ik 2 /deg 
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PARAMETER 

CROSS GEARING 

K AY 

% 

LATERAL ACCELERATION COMPENSATION 
LATERAL ACCELERATION NOTCH FILTER 

LATERAL ACCELERATION 
STRUCTURAL FILTER 

PEDAL GEARING 

ROLL TO YAW COMPENSATION 

YAW RATE STRUCTURAL FILTER 

YAW RATE WASHOUT 

YAW RCV GEARING 

YAW RCV SERVO 


VALUE UNITS 

3.05 DEG/IN 

.83 DEG/FT/SEC 2 

.669 DEG/DEG/ SEC 

(1 + .25S)/ (1 + .125S) 

(S 2 + 2.01062S + 404.26)/ 

(S 2 + 20.1062S + 404.26) 

164/ (S 2 + 4S + 164) 

-7.07 DEG/IN 

.58/(1 + .0435S) 

1/(1 + .0435S) 

S/(l + 3. 11S) 

.35 IN 2 /DEG 

1/(1 + . 0175S) 


FIGURE 2.5-7 

DIRECTIONAL STABILITY AUGMENTATION PARAMETER VALUES 






The propulsion syszm of the YAV- 88 aircraft; consists of the tT4Q2~KI~4G4 
vectored thrust engine, ’fcuAlt: by tolls Boyce Limited, Aero Sivisiap-Btistol, 
ac,d as ait induction system. The engine is similar to the F402-M-402 in the 
M-Ba vtzh a number of design changes to improve performance sad. satisfy 
iacreasec aircraft power extraction requirements. thme changes istiade zero 
degree scarf angle frost aoorles and an Saproved gearbox, •• 


The air induction syszm retains the a?~ 8A basic design concept but 
iacorporates refinements to tasprove ialat oatfotmnce and reduce- distortion. 
This section describes the physical and .operational features of these propul' 
sion system components. 


2*6. 1 YF402~tH«-40d £236 X8S - the TT4O2-RK-404 ehgio* shown in Figure 2.6-1 is 

based on long and extensive development of the fegmm .engine fasti y which - 
includes Pegasus 6, 10 md 11 engines, the engine provides lift thrust for 
takeoff- and landing, cruise thrust for conventional wiagoorae flight, deflected 
thrust for inflight maneuvering and compressor bleed air for the aircraft 
reaction centre I systsa, This versatility, has been achieved through the use 
of a -unique' thtust-vettoring nettle system on. a conventional turfaofaa. aagise 
designed for. l/STQi operation. The mtzle syszm can direct the engine thrust 
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2.6. 1.1 General Description - The YF402-RR-404 engine retains the same basic 
dimensions and features of the F402-RR-402. The engine consists of fan and 
compressor assemblies, an annular combustor, two turbine assemblies, and 
the thrust-vectoring exhaust nozzle system. The significant engine features 
are : 

o Non-mixed-flow turbo fan 

o Two-spool construction - High Pressure (HP) and Low Pressure (LP) 
rotors are contra-rotating to minimize gyroscopic effects. Four 
main bearings are used to support the two rotors (two bearings per 
rotor) . 

o Three-stage fan - The first two stages are overhung from the front 
fan bearing. Inlet guide vanes are not used, thus reducing Foreign 
Object Damage (FOD) susceptibility and eliminating the requirement 
for anti-icing. 

o Eight-stage high-pressure compressor - Incorporates variable inlet 
guide vanes and interstage bleed valves. Aircraft services bleed 
is extracted from the 6th stage, and compressor discharge air (8th 
stage) is provided for the aircraft reaction control system (on 
demand) . 

o Annular combustor with J-tube vaporizing fuel injection 

o Combustor water injection 

o Two-stage air cooled HP turbine 

o Two- st age LP turbine 

These engine features are illustrated in Figure 2.6-2 and significant 
design characteristics are tabulated in Figure 2.6-3. 

An uprated gearbox driven by the high pressure rotor is runted on tor 
of the engine. Drive pads are provided for the engine fuel pump, HP tach- 
ometer generator, inlet-guide vane control, gas-turbine starter, and (2) 
hydraulic pumps. A shaft geared to the low-pressure rotor is used to drive 
the LP tachometer-generator. 

The engine fuel control is a hydromechanical unit mounted to the accessory 
gear box. The fuel control unit incorporates low pressure and high pressure 
fuel pumps, a low-pressure lotor speed governor, and acceleration and flow 
control units. 

The engine oil system is self-contained and is of the scavenge-retum 
type, incorporating a fuel-cooled oil cooler. The oil tank is mounted on 
the engine. 
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FIGURE 2.6*2 

PEGASUS ENGINE DESIGN FEATURES 


Fan stages 

Fan airflow 

Fan pressure ratio 

Fan speed 

Compressor stages 

Compressor airflow 

Compressor speed 

Compressor pressure ratio 
Overall pressure ratio — 

Bypass ratio 

HP turbine stages 

LP turbine stages 


3 

430.7 Ibfeec 

129 

6,500 rpm 
8 

184.2 Ib/sec 
11,000 rpm 
6.0 
13.7 
1.34 
2 
2 
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FIGURE 2.6-3 

YF402-RR-404 SPEC ENGINE CHARACTERISTICS 
100% Corrected Fan Speed 
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The engine is equipped with four rotatable nozzles, two on each side 
of the engine. All four nozzles are a fixed area, convergent type. The 
nozzles are mechanically interconnected and are positioned simultaneously 
by an air actuated drive mechanism which provides the desired degree of 
thrust vectoring. The front (cold) nozzles are airframe mounted on each 
side of the fan delivery duct and are interchangeable between sides. They 
incorporate a zero degree scarf angle. This configuration was developed in 
an extensive test program conducted jointly by Rolls-Royce (RR) and MCA1R. It 
was originally conceived to reduce an adverse exhaust plume/wing/pylons 
interaction that exists with the AV-8A front nozzles. The zero scarf nozzle 
successfully demonstrated a reduction in these adverse interactions in full 
scale tests at NASA Ames Research Center. In developing this configuration 
to a flight worthy nozzle for the YAV-8B, RR has also improved the nozzle 
efficiency and flow characteristics. Test stand calibration of the first 
YAV-8B engine, conducted by Roll -Royce during April 1978, demonstrated a 
200 lb. thrust improvement at the Short Lift Wet rating. 

The rear (hot) nozzles are mounted on the turbine exhaust duct, are 
also interchangeable and are unchanged from the AV-8A/F402-RR-402 . The 
combined effective flow areas are 389 sq. in. for both front nozzles and 
530 sq. in. for both rear nozzles. Each nozzle incorporates two airfoil 
section turning vanes spaced equidistant across the nozzle exit. 

The nozzle control system consists of a nozzle control lever mechanically 
connected to duplicate independent air motors which are, in turn, connected 
to the nozzle drive system. The nozzle system is shown in Figure 2.6-4. 

The YF402-RR-404 retains the current version of the Dowty hydromechanical 
fuel control unit, including all approved Reliability and Maintainability 
modifications . 

Other features incorporated in the YF402-RR-404 include: 

New bulkhead to match the new wing lower moldlines 
Fire extinguisher pipes deleted 
Simplified electrical harness 

External dimensions and details of the engine are presented in Figure 
2.6-5. 


2.6.2 ENGINE OPERATION - Engine operation involves three basic modes; engine 
starting, V/STOL flight, and conventional flight. The engine is controlled 
through two primary systems. These are the engine fuel control and the 
nozzle vectoring control. The fuel control regulates engine speed and 
temperature (and therefore thrust) while the nozzle control system controls 
the direction of the engine gross thrust vector. 
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PEGASUS NOZZLE DRIVE SYSTEM 
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2. 6. 2.1 Fuel Control System - The YF402-RR-404 engine is controlled by a 
hydromechancial fuel control system which governs fan speed and maintains 
operation within selected design limits. The control system components 
are: fan speed governor, limiters, pressure drop regulator for the main 

metering valve, acceleration control unit, air bleed .reset unit (which 
provides additional fuel flow whenever reaction control system (RCS) bleed 
is demanded) , water bypass solenoid, fuel gallery and check valves , igniter- 
jets and main-jets. Engine operating limits are maintained by a set of 
limiters which operate on engine internal pressure and temperature signals. 

The control system also has the capability to override the limiters and to 
select an alternate manual fuel control circuit which serves as a backup 
control system. Additional engine control is provided by the Inlet Guide 
Vane (I GV) control unit, which regulates compressor XGV position and Blow-off 
(bleed) Valve operation. The engine controls are illustrated in Figure 2.6-6. 

The governing function of the control operates above 85Z fan speed by 
controlling the fuel metering orifice area as a function of fan speed, thus 
maintaining a constant relationship between engine speed and power lever 
angle. Acceleration and deceleration circuits override the governor when 
an error exists between the actual fan speed and the power lever commanded 
fan speed. The Acceleration Control Unit (ACU) schedules maximum fuel flow 
as a function of compressor discharge pressure (P 3 ) , and fan delivery pres- 
sure (?13) through positioning of the main metering valve. At less than 85Z 
fan speed the throttle position directly selects a metering orifice area 
independent of fan speed, thus directly controlling fuel flow. In this 
regime, the fuel flow is independent of altitude and ambient temperature. 

The engine speed is determined by the fuel flow with no control or feedback 
into the governing system. Fan speed will vary at a given throttle position 
as ambient conditions vary. 

o Limiters - The limiters sense pressure and temperatures and compare 
the observed values against reference limit values (termed datums).. When 
a limit is reached or exceeded the fuel control reduces fuel flow as 
necessary to maintain engine operation within the specified limits. The 
limits used to control the engine are exhaust gas temperature (EGT) , 
compressor discharge pressure, and HP compressor pressure ratio. 

The EGT limiter is an electronic control unit which senses EGT. When 
the EGT exceeds datum level, the limiter acts through the fan speed governor 
to reduce fan speed until the corresponding EGT is within the proper limits. 
The EGT limiter has three temperature datums. Short lift Wet (SLW), Short 
Lift Dry (SLD) and Maximum Thrust. The specified datum to which the limiter 
operates is dependent on the aircraft mode of operation. The EGT datums are 
interlocked with the nozzle position and landing gear electrical circuits 
so that the V/STOL datums, SLW and SLD, can only be used when either the Q 
landing gear is selected down or the nozzles are rotated to greater than If. 
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ENGINE CONTROL DESCRIPTION 
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The Pressure Hatic Ljmiter (PEL) maintains stable engine operation 
at high altitude/low speed conditions by limiting maximum engine speed, and 
therefore maximum airflow. Corrected fan speed is limited to a nwr-imimi of 
1012 through limitation of the pressure ratio across the HP compressor. 

This is achieved by the PEL, which senses fan and HP delivery pressure, 
bleeding fuel from downstream of the main metering valve back to pump inlet. 

A signal from the air data commit er renders the PEL inoperative at altitudes 
below 10,000 ft where corrected speed limiting is neither needed nor desired. 

l..e Compressor Pressure Limiter (CPL) senses HP compressor exit 
pressure and operates at approximately 255 psig to restrict the P 3 pressure 
signal to the ACU, thus controlling the metered fuel supply. The effect is 
to limit the maximum combustion chamber pressure and also the LP shaft 
torque to values which permit the design cyclic fatigue life to be achieved 
on these items. 

In case of emergency, the limiters, with the exception of CPL, can be 
overriden by the plxot. A switch on the throttle quadrant cuts the pressure 
ratio and exhaust gas temperature limiters out of the control circuit. The 
switch can be actuated either by hand or by pushing the throttle hard against 
the forward throttle stop. Once the limiters have been overridden, it is 
necessary for the pilot to manually control engine speed and temperature 
through the throttle. 

0 Inlet Guide Vane Control - The Inlet Guide Vane (IGV) Control 
regulates one HP compressor variable inlet guide vanes. They are operated 
on a linear schedule against corrected compressor speed from 440° to -2°. 

The IGV control uses fuel as its hydraulic medium but is otherwise completely 
separate from the main fuel system. 

o Blow Off Valve Control - The Blow Off Valve (BOV) Control regulates 
the HP compressor bleed valves oq the 5th stage which bleed approximately 
102 of the core airflow into the plan . chambe These are closed by a 
rotary valve mounted on the IGV operating shaft. The BOV closes/opens at 
18° IGV angle which corresponds to approximately 702 HP RPM/552 fan RPM. 

With the BOV open, additional HP compressor surge margin is obtained at 
low speeds to assist starting and acceleration. 

o Water Injection - Water can be injected into the combustion chamber 
to augment thrust for V/ST0L operation. The control system resets to 
higher speed and temperature limits since the use of water allows increased 
turbine inlet gas temperatures without Increasing metal temperature above the 

corresoondlns dry value. The water flow rate is approximately 350 pounds per 
minute. 

o Manual Fuel Control - In the event of a malfunction in the primary 
fuel control, the pilot can select manual fuel control operation. In this 
mode all governing, acceleration, deceleration, and limiting functions are 
bypassed and fuel flow is directly and mechanically controlled by the power 
lever. All limits must be pilot observed. 
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2. 6. 2. 2 Nozzle Control system - The vectored thrust feature of the Pegasus 
engine enables the pilot to rotate the engine gross thrust vectoi to the 
optimum angle for a given flight condition. Thrust vectoring is accomplished 
through the simultaneous rotation of the four nozzles within the range of 

0° to 98.5° relative to the engine centerline. This rotation is initiated 
by, and is proportional to, the movement of the nozzle control lever. This 
mo - iment is transmitted via pulley and cable to duplicate air motors -/hicb are 
driven by sixth stage compressor bleed air. The air motors drive the nozzles 
through a system of shafts, gears, and chains. The nozzle control lever is 
the only additional cockpit con' 1 required for V/STOL operation; the only 
additional cockpit instrument is che gauge which displays the angular position 
of the nozzles. 

2. 6.2. 3 Operating Limits - The YF402-RR-404 engine control and time limits 
are shown in Figure 2.6-7. The time limitations shown are cumulative and 
include any time spent at higher power settings. The pilot must .nanually 
control the throttle so as not to exceed the fan «peed limits at high 
ambient temperatures while using the Short Lift ratings at zero or low 
bleed rates. The Normal Lift and Maximum Continuous ratings also require 
pilot control over angine speed and temperature. Selection of dry limits 

by the fuel control is dependent upon the operating mode of the aircraft, i.e., 
V/STOL or conventional flight. This information is supplied by electrical 
signals taken from the nozzle and main landing gear positions. VJr t limits 
are actuated by pilot selection of water injection. Below 10,000 ft the 
pilot must manually control engine speed to stay within corrected speed 
limits since the PRL is inoperative. 

For flight test purposes, the YF402-RR-404 engines have be^n cleared to 
operate at Combat EGT levels for wingborne flight and for extended periods 
within Normal L? ft Dry for hover. Combat EGT levels of 650°C or 675°C are 
available for a maximum of 2.5 minutes in wingborne flight. This is 
achieved through the use of special EGT Limiter datum temperature selector 


Rating 

Fan speed 

limit 

(% N F ) 

Corrected 
fan speed 
limit (% N f / 

Maximum 
exhaust gas 
temperature 
limit 

rc> 

Compressor 

discharge 

pressure 

limit 

(psig) 

Cumulative 

time 

limir 

Short lift (wet) 

107.0 

106.5 

745 

255 

15 sec 

Short lift 

103.5 



715 


15 sec 

Normal iift (wet) 

104.5 



720 


1-1/2 mm 

Normal lift - * 

100.0 



695 


2-1/2 min 

Maximum thru 

, 95.5 

■ 

101.0 
•tx»e 
10.000 ft 

610 


15 min 

Maximum cont 

89.0 

■ 

540 


None 


‘Includes time at all higher ratings 
“May also be used in flight as a combat power setting 


FIGURE 2.6*7 

YF402-RR-404 ENGINE OPERATING LIMITS 
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plugs, which substitute either 650°C or 675®C for the normal Max Thrust 
limit of 600 e C. When either of the special plugs is installed, the Max Thrust 
EGT limitation must be manually controlled by the pilot. In addition, when 
operating in the Combat EGT regime, uhe fan speed must be limited to 1001 
at airspeeds between 250 to 450 knots, and to 95. 5X at airspeeds above 
450 toots. An engine life penalty is also assessed ' tech flight using 
Combat EGT. 

For the purposes of achieving steady state conditions during sustained 
performance hover activity, the YF402-RR-404 has been cleared to operate for 
4 minutes at exhaust gas temperatres below 685 C, in lieu of the normal 2.5 
minute limitation within Normal Lift Dry (NLD). Four minutes at 685*C 
’■ suits in the same number of Engine Life Recorder counts as does 2.5 minutes 
the NLD EGT limit of 695 °C, therefore no engine life penalty is associated 
rfith this extension. 


2. 6.2. 4 Engine Starting - The YAV-8B propulsion system is completely self- 
sufficient for engine starting. An engine mounted gas turbine starter 
(GTS) is started using internal battery power. The TAV-8B incorporates an 
Improved Lucas GT C -. which is derived from the Lucas Mx I and Mk II GTS units 
used on the AV-8A. The aircraft boost pumps are actuated by switches on 
the righthand instrument panel. Opening the LP fuel valve then allows 

fuel to be delivered to the engine LP fuel pump. Moving the START switch 
to ON opens a solenoid valve to deliver fuel to the GTS and arms the starting 
circuits. The momentary movement to the START position initiates the GTS 
and engine starting cycle. 

The GTS is started by a 3.0 horsepower electrical motor powered by an 
aircraft battery. The electric motor automatically disengages when the 
CTS reaches self-sustaining speed. For engine starting, the GTS power turbine 
drives the engine through the auxiliary gearbox. Rotation of the HP spool 
allows fuel to be delivered to fuel-actuated engine accessories such as 
the IGV control. After the GTS has started, the HP fuel valve must be 
opened by advancing the engine throttle to the IDLE position. This allows 
fuel to be delivered through the fuel control units to the engine combustor. 

By the time the GTS output shaft *nd the HP rotor of the engine have acceler- 
ated to 3200 RPM, the engine is -f- sustaining and starter assist is no longei 

required. The GTS is then automatically shut down by shutting off its 
fuel supply. The engine continues to accelerate up to its idle speed of 
approximately 251 fan speed. 

tn order to shut dowu the engine, the throttle is simply retarded to 
the cutoff position, which closes the HP fuel valve and stops the fuel flow. 

2. 6. 2. 5 V'STOL Operation - In the V/ST0L flight mode, the EGT limiter is 
reset from the wingborne flight limit to the lift rating limit. The thrust- 
vectoring exhaust nozzles are down and compressor bleed air is provided to 
the aircraft reaction control system. 
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Whenever the nozzles are deflected more than 16° below the aft position 
or the main landing gear is locked down, the exhaust gas temperature limiter 
is reset to the Short Lift Dry value to allow operation at higher engine 
speed. Thrust augmentation is available by water injection into the combus- 
tor if the fan speed is at least 95%. Selection of water injection by the 
pilot turns on the water pump and resets the fan speed governor to a higher 
limit. The EGT limber is reset to a higher value as soon as the water 
begins to flow. Automatic reversion to the Short Lift Dry rating will occur 
if the water supply is exhausted, the fan speed drops below 75%, or the pilot 
shuts off the water injection system switch. 

Automatic engine control within the applicable speed and temperature 
limits is provided only in the Short Lift ratings; the pilot must manually 
control the engine tc the Normal Lift speed and temperature limits. 

Rotation of the nozzles beyond 36° fully opens a butterfly valve in the 
aircraft reaction control system duct that makes .compressor discharge bleed 
air available for attitude control during takeoff and landing. Maximum air* 
flow is available for nozzle deflections beyond 15 . When the lift ratings 
are used, total thrust (including the thrust from the reaction control system) 
is maintained approximately constant over a range of bleed rates. (The 
extent of this range is dependent upon ambient temperature and power setting) . 
This is accomplished by increasing fuel flow and the turbine inlet tempera- 
ture as the bleed rate is increased. An aix-bleed-reset within the fuel 
control unit provides the additional fuel flow during acceleration whenever 
bleed air is being extracted. As bleed is increased, the exhaust gas tempera- 
ture will eventually reach a selected limit. Further increases in bleed then 
cause a reduction in thrust at the limiting temperature. When the nozzles 
are rotated aft, reaction control bleed is shut off. Also, when the nozzles 
are rotated aft to their conventional flight positions and the main landing 
gear is retracted, the EG1 limiter resets to conventional flight limits and 
the engine is controlled to the Maximum Thrust rating. 

2. 6.2. 6 Conventional Flight Operatior - Engine operation in the conventional 
flight mode is similar to that of other engines. The highest power setting 
in normal operation is a 15 minute Maximum Thrust rating. The YF402-RR-404 
engine-* have been cleared to a 2.5 minute Cor ^at rating for flight test. 

Thrust vectoring may be used as an aid to aircraft maneuvering, within th rt 
flight envelope restrictions. The nozzle system can be rotated through 98.5 
to obtain reverse thrust either in flight or on the ground. 

2.6.3 AIR INDUCTION SYSTEM - The air induction system for the YAV-8B is the 
direct result of an inlet development program which follows a minimum change 
philosophy. The basic AV-8A inlet concept has been retained with modifica- 
tions incorporated to achieve desired performance and distortion improvements. 
The inlet development program was directed toward an increase in V/STOL total 
pressure recovery and reduction in high speed distortion at high positive and 
negative angles of attack. Extensive inlet model testing at static, low 
speed and high speed conditions, combined with full scale test results at 
NASA Ames Research Center facilities, has shown that these objectives have 
been achieved with the YAV-8B inlet. 
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2. 6. 3.1 Physical Description - The YAV-8B Inlet is a subsonic bifurcated pitot 
type and incorporates two side inlet ducts each of approximately 190 ° circular 
arc at the highlight plane. The ducts merge two inches ahead of the engine 
face. The main lip incorporates a large contraction ratio (highlight area 
divided by lip throat area, 1.24) where the lip throat is defined as the 
junction between a 2:1 quarter-ellipse lip and the internal duct lines. 

This lip shape is being used on the YAV-8B inlet because of its excellent 
V/STOL and high speed performance characteristics . 

A large auxiliary inlet around the periphery of each cowl provides 
additional inlet area for takeoff and low speed operation. The double- 
door auxiliary inlet configuration is shown in Figure 2.6-8. Each auxiliary 
inlet contains seven passages which converge to form a continuous slot 
before entering the main inlet duct. The passages are separated by fore- 
and-aft struts which connect the front and rear sections of the inlet cowl. 

Each auxiliary inlet passage incorporates two free floating blow- in- doors 
hinged at their leading edges. The door external surface is flush with 
the aircraft moldline in the closed position. 

The boundary layer bleed system, shown in Figure 2.6-9, is similar to 
that of the AV-8A. A flush bleed slot is located on the inboard wall of 
each inlet, immediately aft of the entry plane. The slot is covered by 
a two piece door (split into upper and lower halves; which is hinged at 
its leading edge and Is spring-loaded closed. Boundary layer air from each 
slot is ducted upward to an exit at the rear of the canopy, from which it 
discharges overboard at an angle of approximately 45°. 

The inlet plan and side views previously shown in Figure 2.6-8 illustrate 
the short highly curved nature of the inlet ducts. The duct length is 
roughly equal to the engine face diameter due to constraints set by the 
overall airframe layout. Figure 2.6-10 presents the inlet design parameters. 

2. 6. 3. 2 Operation - The YAV-8B inlet has two separate and distinct operating 
modes, i. e. , V/STOL and conventional flight operation. These modes are 
determined by the positions of the auxiliary inlets and the boundary layer 
bleed doors. Scctic pressure differential across these components 
determines their positions and hence the operating mode. The auxiliary inlet 
doors are fr~ --floating, their positions determined only by static pressure 
different- 1 . The spring-loaded boundary layer bleed doors are positioned 

by a st pressure differential and an opposing spring force. 

o V/STOL Operation - The engine airflow required for vertical or short 
takeoff is supplied at high total pressure recovery through use of a 
large inlet contraction ratio and auxiliary inlets. The auxiliary inlets 
increase available flaw area, thus maintaining low velocity past the main 
inlet lips and minimizing total pressure loss. During V/STOL operation the 
free-floating auxiliary inlet doors open due to the pressure differential 
caused by the duct static pressure being less than ambient pressure. The 
combined effects of spring force and pressure differential hold the boundary 
layer bleed doors closed. Figure 2.6-11 illustrates inlet operation in the 
V/STOL regime. 
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FIGURE 2.M 

GENERAL INLET ARRANGEMENT 



BOUNDARY LAYER BLEED SYSTEM 
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YAV48 

Minimum duct area (ft 2 ) 

Up shape 

Highlight area (ft 2 ) 

Mean lip contraction ratio 
Auxiliary throat area (ft 2 ) 

Auxiliary entry area (ft 2 ) 

Auxiliary contraction ratio 
Maximum cowl area ,‘t 2 ) 

9.7 

2:1 ellipse 
1253 
1.24 
8.41 
1258 
1.49 
17.78 




RGURE 2.6-10 

INLET DESIGN PARAMETERS 



FIGURE 2.6-11 
INLET OPERATING MODES 
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o Conventional Flight Operation - During conventional flight (cruise 
or high speed) the auxiliary inlet doors are closed as a result of the inlet 
duct pressure being higher than ambient* The boundary layer bleed doors 
are positioned as a function of the duct static pressure which is established 
by flight Mach number and engine power setting, the bleed exit static pressure 
(usually ambient), and the opposing spring force. At cruise power settings 
the reduced inlet mass flow ratio results in reduced inlet velocities with 
resultant increased duct static pressure. As the spring force is overcome 
by the pressure differential across the door(duct pressure minus bleed exit 
pressure), the doors open to remove the boundary layer from the inlet flow. 

The boundary layer bleed doors are free to operate individually in response 
to static pressure gradients across the doors during operation at angles of 
attack. The split upper and lower door configuration developed for the AV-8A 
and retained for the YAV-8B is used to prevent recirculation of boundary 
layer air from a high pressure region to a law pressure region, which would 
cause separation and increased inlet pressure distortion. Figure 2.6-11 
illustrates conventional flight operation. 
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The nonlinear program that has been developed to simulate the YAV-8B 
aircraft is named YAV8B. The program computes six degree-of-freedom aircraft 
motion plus a large number of aircraft performance parameters. The program is 
written in FORTRAN in a manner to make it easy to use. Program organization, 
structure and data card input format are described within this section. 

3.1 PROGRAM STRUCTURE 

The YAV8B program is structured into two primary categories: the 

executive routines and the aircraft math model routines. For simplicity, all 
computations are performed at the same iteration time. The program computes 
aircraft motion as well as performance parameters. A great deal of flexibility 
is available to the user with respect to which parameters are to be recorded. 

Figure 3.1-1 diagrams the structure of the YAV8B program. The executive 
(exec) routines include I5O0NCE, 150M5, AC07, CARDS and RTPDATA. The I5O0NCE and 
150MS routines are top level exec routines for the program. The AC07 routine 
is the aircraft math model exec. The CARDS routine reads the input data from 
the EXECUTE card deck and the RPTDATA routine is the exec for the output data 
logic. The program computes in three different modes: INITIAL (first pass 

initialization) , I RESET ( 400 passes for dynamic stabilization) , and OPERATE 
(real time computation) . The run time for each particular case can be 
designated. The iteration time (DT) is defaulted to 50 milliseconds but can be 
changed by the user. 

Initial condition data, as well as the instructions for Cata printout, 
are input to the program using the EXECUTE card deck. The EXECUTE card deck 
loads the appropriate software files and executes the program. Dynamic 
inputs to the aircraft controls, aircraft motion, etc. can be input to the 
program with the I50MS card deck. The I50MS deck modifies the I53MS sub- 
routine. Because the I50MS subroutine is a top-level exec routine, this is 
a convenient location to math model a forcing function. A more detailed 
description of user procedures is provided in Section 3.4. A sample case is 
also provided in Section 3.4. 

The aircraft math model flow diagram is shown in Figure 3.1-2. Each 
subroutine in the aircraft math model is described in Section 3.3 along with 
top-level flow charts. All of the subroutines for the aircraft math model 
can be modified and changed by* the user if required. 

3.2 MATH MODEL EQUATIONS 

The YAV8B program computes aircraft velocities, accelerations and 
position at the aircraft c.g. The airframe is considered to be a rigid body. 

The atmosphere data base is the U.S. Standard Atmosphere, 1962 - sea level to 

100.000 feet. The following paragraphs describe some of the theory and math 
model equations comprising the aircraft simulation model. 

3.2.1 N UMERICAL METHODS - There are three general numeric methods used in 

the YAV8B program that require some discussion. They are integration 
techniques, digital filter techniques, and data table look-up methods. Each 
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FIGURE 3.1-1 

YAV8B PROGRAM STRUCTURE 
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AIRCRAFT MATH MODEL FLOW DIAGRAM 
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n«thod will be addressed briefly in order to familiarize the program user 
with the techniques incorporated into TAV8B. 


3.2. 1.1 Integration Techniques - There are four types of integration schemes 
used in the simulation program: Euler (rectangular) , T us tin, Reddy Parabolic 

and Adams third order. These methods fall into two general classes. The 
first class is the predictor integration scheme. In this class, the most 
currently available input is valid for the beginning of the integration inter- 
val and the output is "predicted" for the end of the Interval. The second is 
the corrector integration scheme. This scheme is used when the input is 
available at the same time that the output is desired. The Adams third order 
method is used primarily in computations for the equations of motion. The 
general form of these numerical integration methods are as follows: 

o Euler (Rectangular) 

X n - x x + X Q _ 1 * At, At - Sample period Zero Order Predictor 


o Tustin 


X_ 


X . + (.5 i + .5 X ) * At 
n-1 n n-i 


o Adams Third Order 

V Vi + <§ *-1 - I V: + u W * “ 


Linear Corrector 


Third Order Predictor 


o 


Reddy Parabolic 

*n ■ Vl + <§ *n + 3 Vl 


12 


K-J * “ 


Parabolic Corrector 


3. 2. 1.2 Digital Filtering Techniques - Digital filters are needed to obtain 
an acceptable model of hydraulic actuators in the control system model and 
for other lead and lag filters used in the simulation model. Two methods are 
used in the YAV8B program: Tustin ' s and Reddy's. The Tustin method is quite 
good for modeling digital filters, as it has a simple relationship between the 
Laplace transform and the Z transform operators ? and Z: 


1 . At (Z + 1) 
S 2 (Z - 1) 


Z - Z Transform operator 


The equivalent computer difference equation for C/R 


1/S would be: 


C ■ C , + [R + R , ] Where R is the present value of the input, 
n n-1 2 n n-l J n r 

Using the same method, the difference equation for the first order filter, 
C/R - 1 /tS+I, can be derived as: 


C 

n 


(1 + ~) 
v At' 


r- <R + R r ) 

n n- 1 


(i - ~) 

2^7 C n-1 » 


(1 *S? 


r - system tl~e constant 


This difference equation would give the output sequence (C n ) of a first order 
lag from the input sequence (Rq) . Note that the coefficient terms involve only 
the constants At and t, and so are only computed once at the beginning of a run. 
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For better accuracy and stability, some filters with short time constants 
are modeled using Reddy's method. For a first order transfer function and 
linear input, the output at the eno of each step is a constant times the 
output at the beginning of each step plus other constants times values 
of the input at the beginning and erd of the step. A first order system 
response to a linear input can be represented by the following recursion 
formula: 


C n * C n-1 (e ‘ At/T > + \tl* (t/ At) (1 - e~ At/T )] 

+ R n _ 1 [t/At - (e' At/T )(l + t/ At) ] 

3.2. 1.3 Data Table Look-Up Techniques - The method used in the table look-up 
program, TBLKP, for even increments is simple and straightforward. It 
requires the program to compute an index and a ratio. The index represents 
the position in the table of dependent variables such that the corresponding 
value of the independent variable is just below (or equal to) the input 
variable. The input variable is bound by this value and the next value of 
the independent variable. The ratio is the percentage distance of the input 
variable between these two values. 

For a single dimension table, the lookup function performs as shown : 

RI -((Input - low bound) (number of points - 1) / (high bound - low bound) ) +1 
Index * integer part of RI 
Ratio * fractional part of RI 

Once the index and ratio are found, it is only necessary to look up the 
corresponding dependent variables and linearly interpolate: 

Output * Table (Index) + [Table (Index + 1) - Table (Index)] x Ratio 

The lookup method for a table with uneven breakpoints i$ similar, but 
each set of breakpoints must be checked to determine wnlcb ones bracket the 
input value. 

Since the values of the independent variable will cover the allowable 
boundaries, the table lookup program does not extrapolate outside the 
boundaries of the tables. If an independent variable exceeds the boundary, 
the corresponding dependent variable boundary value is used in the lookup. 
Finally, to save additional time with either method, the index into the table 
and the ratio between points is saved, so that if more than one table uses 
the same breakpoints, the index and ratio need not be recomputed. 

3.2.2 EQUATIONS OF MOTION - The parameters that determine the aircraft motion 
are described by a set ol six second-order differential equations, three to 
describe the translational degrees of freedom and three to describe the 
rotational degrees of freedom. This system of equations is used to describe 
the aircraft state relative to a fixed coordinate reference frame. 
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The earth is modeled as a spherical surface, rotating relative to a fixed 
inertial frame. Rotation is about the polar axis with no "wobble". The 
acceleration due to gravity (g) is assumed to be constant and always acting 
perpendicular to the ground. 

Four coordinate reference frames are used in describing the aircraft 
position, velocity and acceleration: 

(1) Inertial Axes 

(2) North-East-Down (NED) , or local vertical 

(3) Body Axes 

(4) Stability Axes 

The inertial coordinate frame is an orthogonal triad with its origin at 
the center of earth (see Figure 3.2-1). The first and second components (or 
x and y axes) define the equatorial plane and the z axis is positive toward 
the North Pole. The x axis orientation is initialized at the beginning of 
the mission to pass through the zero . latitude and longitude point. The y 
axis passes through the point of latitude • 0°, longitude “ +90°. After the 
mission begins, the rotation of the earth will cause an earth fixed longitude 
point to Increase in inertia j. longitude (relative to the "fixed in space" 
inertial axes reference frame) . It is in this axis system that accelerations 
and velocities are integrated to obtain aircraft position. 


North Pole 



FIGURE 3.2-1 

INERTIAL COORDINATE FRAME 
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The NED frame is centered at the aircraft c.g. The x axis is directed 
toward true North, the y axis is directed toward true East, and the z ax±E *s 
directed down perpendicular to the ground. Figure 3.2-2 illustrates this 
system. Aircraft translational velocities are transformed from inertial axis 
into this frane for the addition of wind and computation of ground speed and 
flight path angles. The aircraft Euler angles (pitch, roll, yaw) are defined 
from the body axes orientation relative to the NED reference frame. 



D GP79-0269-343 


FIGURE 3.2-2 

RELATIONSHIP OF BODY TO NED AXES 


The equations of motion and their computational arrangement are 
illustrated in Figures 3.2-3 and 3.2-4. 

The translational equations are a summation of the total aerodynamic, 
engine thrust, and RCS forces in the body axib coordinate system. Newton's 
Second Law of Motion is applied to compute the body axis accelerations and 
these are transformed to the inertial frame through a direction cosine 
matrix, as shown on Figure 3.2-3. The acceleration due to gravity is added 
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and the acceleration components are integrated to produce inertial velocity 
and position components. These are transformed back to the NED axes for the 
introduction of wind components and to the body axes for the computation of 
the body axis velocity components and the velocity vector angles (a, 8, y, c) . 

The rotational equations of motion are derived by summing the moments due 
to aerodynamics, engine thrust, and ECS effects. The body axis rotational 
accelerations are computed from these moments, the aircraft moments of 
inertia, and the cros- coupling terms. The body rates are determined by 
integrating the rotational accelerations. 

The attitude angles of the simulated aircraft are calculated using the 
quaternion method. The quaternion method uses four parameters for specifying 
the orientation of a coordinate system. The quaternion rates are calculated 
from the body axis angular rates and then integrated to update the quater- 
nions. The earth- 1 o-b ody axis transformation matrix is computed directly 
from the quaternions, thereby defining the aircraft orientation with respect 
to the fixed frame. A correction factor is applied to compensate for inte- 
grator drift and to ensure chat the coordinate system remains orthogonal. 

This correction factor is based upon the constraint that for an orthogonal 
system, the sum of the squares of the quaternions must be unity. The quater- 
nion method offers the advantage of all owing an all-attitude orientation that 
is free from discontinuities. 

3.2.3 AERODYNAMIC EQUATIONS - The aerodynamic forces and moments acting on 
the aircraft are fully described by algorithms for three longitudinal and 
three lateral -directional coefficients. Ths longitudinal terms are the 
coefficient of lift (Cl) or normal force (Cg) , drag (Cq) or axial force (C A ), 
and pitching moment (Cg) . The lateral-directional terms are the coefficients 
of side force, (C^) , rolling moment (Ci) , and yawing moment (C n ). 

3. 2. 3.1 Aerodynamic Reference Systems - The aerodynamic data base is divided 
into two separat' and distinct regions: low speed (M ^ 0. 3) and high speed 
(M > 0.5). The data for the low speed region are derived in the body asds 
system. The aircraft body coordinate system is shown in Figure 3.2-5a. It 
has its origin at the aircraft's center of gravity. The x axis is directed 
toward the nose of the aircraft coincident with the Fuselage Reference Line 
(FRL) . The y axis is directed toward the right wing of the aircraft (ortho- 
gonal to the x axis) . The z axis is directed downward, perpendicular to 
the X£ - y B plane. 

The data for the high speed region are referenced to the stability axis 
frame. Figure 3.2-5b illustrates the stability axis system. It has its 
origin at the point of reference of the aerodynamic data. The x axis is 
formed at the projection of the body axis on the velocity vector (rotated 
through the aircraft angle of attack a ) . The y axis is directed out the 
right wing of the aircraft and is oriented perpendicular to the aircraft 
plane of symmetry. The z axis is directed downward, perpendic^.ar tc the 
x-y plane. 
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For ground hover and low speed flight, another reference frame is 
used. This reference frame is a blehd between the body axis and stability 
axis reference frame. The ground effects are introduced into the simulation 
for altitudes of approximately 7C feet or less. The longitudinal reference 
attitude for ground effects is as follows: 

*'« a * + e * K 2 for 0.05 < Veq < 0.10 

Where » (Veq - .051/20. and Kj « (1 - 

for Veq < 0.05, a' » 6 and for Veq > 0.10, a' * a 

The FORTRAN name for a' is THAI..P . 

The parameter Veq is defined as the square root of the ratio of the freestream 
dynamic pressure to jet dynamic pressure: 



The lateral -directional reference angle is defined as below: 
8* • Tan -1 (v'/u’) 

Where u' * u/cos a and ▼' * v * cos $ 

The FORTRAN name for 8' is BETAP. 
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3. 2. 3. 2 Aerodynamic Algorithms - Each component of the respective coefficient 
is summarized for the low speed and high speed data in the following equations. 
The FORTRAN name of each component is also included in parenthesis. 

o Low Speed Aero (M ^ 0.3) - Body &xis System 

o Normal Force 

h ■ [ S + '■ > + i V 1T 

0 o„ o_ FJX 


H 


+ (AC, 


N. 




) 

\ 


+ C 


Nx 


* 6 


GEAR 


+ C N 


* 6 


GEAR 


6 LID 


LID 


+(C N * < + S. * d) * * + AC N 1 * i * S + 

q a T POWER GE 


Where: 


- Baseline normal force coefficient as a function of 
a, (CNBASE) . 


AC 


N, 


H 


- Increment or normal force coefficient due to stabilator 
deflection as a function of a_ TT , 0 T , Veq and altitude, 
(DCNSTAB) . XAiL J 


AC. 


N. 


- Increment of normal force coefficient due to left and 
_ right flap deflection as a function of a , 0 T and Veq, 
L,R (CNFLAP). J 


AC, 


N. 


- Increment of normal force coefficient due to flap-jet 
FJI impingement as a function of <5 , 0 , a and Veq, (CNFJI). 

F J 


AC - increment of normal force coefficient due to left and 

right aileron deflection as a function of a, (DCNAIL) . 


N 


" S V36GEAR’ Where 6 GEAR is the normallzed 8 ear 
GEAR deflection, (DCNGEAR). 


LID 


- SC.,, , where 6 TTT( is the normalized LIDS deflection, 

N/ 36LID IL 

(DCNLID) . 


AC N - Increment of normal force due to ground effects as a 
GE function of altitude, a, 0j and engine thrust, (CNGE). 

norma ^ f° rce coefficient due to pitch rate, CCNQ) . 


N. 

a 


- normal force coefficient due to angle of attack 

rate, ( CNALPD) . 


AC W - ^-icrement of normal force coefficient due to power effects 
as a function of o, 0j and Veq, (DCNPOW) . 


POWER 
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C - Baseline axial force coefficient as a function of 

o o, tCABASE) . 


AC - Increm en t of axial coefficient due to stabilator 

5 h deflection, CDCASTAB) . 

AC^ - Increment of axial coefficient due to left and right flap 

d deflection as a function of a and Veq, (CAFLAP). 

F L,R 

AC - Increment of axial force coefficient due to flap-j er 

rJI impingement as a function of 6^, 6^, and a, (CAFJI). 

AC. - Increment of axial force coefficient due to left and 

o. right aileron deflection as a function of a, 'DCAAIL) . 

C^ - 3C A /&6 r; irAB where d ^ ^ is the normalized gear aeflec- 

5 GEAR tion, ''DCAGEAR) . 

C. - 3C./36. __ where 6 T __ is the normalized LIDS deflection, 

A, A LID LID 

LID (DCALID) . 

AC - Increment of axial force due to ground effects as a 

GE function of a, Qj, altitude and engine thrust, ^CAGE). 

AC - Increment of axial force coefficient due to power effects 

“POWER as a function of 6j, Veq and a, (DCAPOW) . 
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o Pitching Moment 
PM 


0F fvo.i QiAury 


[ ^ + uC 

ju m. 

o c 


+ (AC 


H 


nu 

c. 


+ AC ) 4* LC 


xn f 


“fji 


(AC 




+ AC 






) + C 


m. 


*R 
* 6 


GEAR + C m t 


* 6 


LID 




GEAR 


LID 


(C 


m 


q + C * d) * -==- + AC 

2V„ 


m 


“power 


] * q * S * o +■ AC 


AC 


GE 


“TP 


Where: 


m 


AC 


AC 


m. 


L,R 


AC 


“FJI 


AC 


m, 


Vr 


m. 


GEAR 


m. 


LID 


- Baseline pitching moment coefficient as a function of 
a, (CMBASE) . 

- Increment of pitching moment coefficient due to stabilator 
deflection as a function of a TAIL , Q j> Veq and altitude, (DQIS 1 

- Increment of pitching moment coefficient due to left and right 
flap deflection as a function of o and Veq, (CMFLAP). 

- Increment of pitching moment coefficient due to flap-jet 

impingement as a function of a and Veq, (CMFJI) . 

- Inc. dent tf notching moment coefficient due to left and 

right .t reflection as a function of angle of attack, 

(DCMAII) . 

- where 6^^ is the normalized gear deflec- 


tion, (DCMGEAR) . 

W m /,< 1IDS Vl “ r ' 
tion, (DCMLID) . 


- dC /36,___ where 6 TTT . is the normalized LIDS deflec- 
XU LXUo LID 


AC, 


M 


GE 


- Increment of pitching moment due to ground effects as 
a function of o, altitude and engine thrust, (CMGE) . 


C 

m 

q 




- 3C /3q, pitching moment coefficient due to pitch rate, (CMQ) . 

m a 

- SC /da, pitching moment coefficient due to angle of attackj (CM/ 


AC 

“power 


- Increment of 
effects as a 


pitching moment coefficient due to power 
function of a, Veq and 6 , (DCMPOW) . 

J 


AC 

“FP • 


- Increment of pitching moment due to flat plate drag 
motion for u < 50 KTAS. Simulates aerodynamic damping 
in hover and low speeds, (TEMPX3) . 


3-15 



ORIGINAL PAGE IS 
OF POOR QUALITY 


MDC A7910 
Volume I 


o Side Force 


SF - [C_ + (Ac + Ac ) + C. 

o *6 6 6, 

*1 ' 


* 0„ + ACy 


R 


POWER 


] * q *S 


AC + ACL 
GE PHI 


+ 


OGE 


Where: 


V 


- Baseline side force coefficient. For u > 50 KTAS 
defined as C * B, side force due to sideslip angle 
Y S 

as a function of a, B, and average flap deflection, 
(CYBASE). For u < 30 KTAS, C^ is 0. 


^R 


AC^ - Increment of side force coefficient due to left and 

’ right aileron deflection as a function of a, (DCYAIL) . 

- side force coefficient due to rudder reflection, 

'R (CYDR>. 

- Increment of side force coefficient due to power effects 
POWER as a function of Veq,a and B for u > 50 KTAS, (CYPOWER). 

- Increment of side force due to ground effects as a 

GE function of a, altitude, engine thrust, 6j and 6 (SFIGE) . 

- Increment of side force due to roll attitude In 

PHI ground proximity (h < 70 ft.) as a function of altitude, Veq, 

d, 6', u, © and engine thrust, (YLSFTP) . 


s. 

G 

ac y 

ac y 


ac y. 


AC - Increment of side force as a function of a, 0 , engine 

Y 0GE thrust, and Veq for u < 30 KTAS, (SFOGE). For u > 50 

KTAS, ACL is 0. 

X OGE 


o Rolling Moment 

EM - [C, + (AC. + AC, 

~o 6 a 6 

+ AC, + (C, * r + C 

*TJI r 


) + C. 


* S R + LC l. 


+ AC. 


J R 


* n) * —— 

0 P ' 2V 
p- V T 


POWER 


] * q * S * b + AC 


PHI 


+ AC, + AC, + AC- 

FF ^GE OGE 
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Baseline rolling moment coefficient. For u > 50 KTAS 
defined as * 8, rolling moment due to sideslip 

angle as a function of a, 6 and average flap deflection, 
(CLLBASE). For u < 30 ITT AS, is 0. 

o 


AC, 


a l,r 


C 


£ 


Increment of rolling moment coefficient due to left and 
right aileron deflection as a function of a, (DCLLAIL) . 

30^/36^, rolling moment coefficient due to rudder 
deflection, (CLLDR) . 


AC 


£ power 


Increment of rolling moment coefficient due to power effects 
as a function of Veq, a and 6 for u > 50 KTAS, (CLLPOWR) . 


AC 


£ 


Increment of rolling moment coefficient due to asymmetric 
flaps, (DCLLF) • 


AC 


FJI 


Increment of rolling moment coefficient due to flap -jet 
impingement , (DCLLFJ1) . 


C 


£ 


r 


C 


3C^/3r, rolling moment coefficient due to yaw rate as a 
function of a, (CLLR) . 

3C./. , rolling moment coefficient due to roll rate as a 

X, op 

function of a, (CLLP). 


AC 


GE 


Increment of rolling moment due to ground 
a function of a, altitude, engine thrust, 
(KMIGE) . 


effects as 
Veq, 6 j and B ' , 



AC 


£ 


FP 


Increment of rolling moment due to roll attitude in 
ground proximity as a function of altitude, Veq, 4, fjl, 
u, 9 j and engine thrust, (DRMPHI) . 

Increment of rolling moment due to flat plate drag 
motion for u < 50 KTAS. Simulates aerodynamic damping 
in hover and low speeds, (TEMFX5) . 



Increment of roiling moment as a function of 8', a, 
Veq, 6 and engine thrust, (RMDGE) . For u > 50 KTAS, 

AC J - 0. 

^OGE 
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n x 

0 


- n. R 
°R 


+ AC + AC 


“power n 6. 


+ (C * r + C * p) * rsH *q*S*b + AC + AC + AC + C 


C - Baseline yawing moment coefficient. Por u > 50 KTAS 

n o defined as C. * 6, yawing moment due to sideslip angle 

8 

as a function of a, 8 and average flap deflection, 
(CLNBASE). For u < 30 KTAS, C is 0. 

“o 


Increment of yawing moment coefficient due to left and 
right aileron deflection as a function of a, (DCLNAIL) . 


- 3C /3<5 , yawing moment coefficient due to rudder 
i. n R 

°R deflection, (CLNDR1 . 

Z - Increment of yawing moment coefficient due to power 

“POWER effects as a function of Veq, a and 8 for u > 50 KTAS, (CUTPOWR) 


Increment of yawing moment fioeff ic lent due to asymmetric 
flaps. (DCLNF). 

Increment of yawing moment coefficient due to flap-jet 
impingement, (DCLNFjl) . 

yawing moment coefficient due to yaw rate as a 
function of a, (CLNR) . 

3C /- , yawing moment coefficient due to roll rate as a 
n op 

function of a, CCLNP) . 

•Increment of yawing moment due to ground effects as a function 
altitude, engine thrust, Veq, 0j and 6' , (YMIGE) . 

Increment of yawing moment due to roll attitude in giound 
proximity (alitltude < 70') as a function of altitude, 

Veq, 41 , 8', u, 9^ and engine thrust, (DYMPHI) . 

Increment of yawing moment due to flat plate irag 
motion for u < 50 KTAS. Simulates aerodynamic damping 
in hover and low speeds, (TEHPX6) . 

Increment of yawing moment as a function of S’, a, Veq, 
and engine thrust, (7M0GE) . For u > 50 KTAS, 

AC « 0 
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o High Speed Aero (M .> 0.5) - Stability Axis System 
o Lift 

F. = [C + (ACL + AC. ) + (AC t + AC ) + AC ] * q * S 

L o L 5_ L 6_ 6 L 6, ^POWER 

F 1 F R *1 *R 


+ ac t 


'STAB 


Where: 

C L 

o 


- Baseline lift force coefficient as a function of a 
and M, CCLBASE) . 


AC^ - Increment of lift force coefficient due to left and 

L & right flap deflection as a function of a and M, 

f L,R (DCLFL and DCLFR) . 


AC T - Increment of lift force coefficient due to left and 

right aileron deflection as a function of M, (DCLAL 
L 9 R and DCLAil) • 


A., 


AC, 


^OWER 

ac l 

STAB 
o Drag 


- Increment of lift force coefficient to power effects 
as a function of 6^, a and Veq, (DCLPOW) . 

- Increment of lift force due to stabilator deflection, 
as a function of a and M, (CLSTAB) . 


F D • + ( 4C I,, + iC D ) + <4 S, + * + 4C I 


F 

L F R 

+ C D *i 6 rI + aC d +C D 

% \ear 


*L *R 


STAB 


* 6 


GEAR + 


* 6 lid ] * * * S 


LID 


Where: 


- Baseline drag coefficient as a function of a, M and 
C , (CDBASE) . 

L o 


AC, 


L,R 


AC, 


D. 

o , 


^.R 


- Increment of drag force due coefficient due to left and right 
flap deflection as a function of a and li, (DCDFL and DCDFR) . 

- Increment of drag force coefficient due to left and right 
aileron deflection as a function of a and M, (DCDAL and DCDAR) . 
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STAB 


- Increment of drag force coefficient due to stabilator 
deflection as a function of a and M, (CDSTAB). 


R 


- drag force coefficient due to rudder deflection 

as a function of a and M, (CDDR) . 


AC. - Increment of drag force coefficient due to pylons and 

STOW gun pods, (HCDSTOW) . 


GEAR 


- * w h®re 6 / ^ rA1> is the normalized gear deflec- 

tion. SEAR 


C_ - 3C_/.. , where 6. _ is the normalized LIDS deflection, 

o Pitching Moment 

PM- [C + (AC + AC ) + (AC + AC ) 


k R 


U U 


+ AC 


“stab 


(C 


q + C * o) * ~ + AC 
™* 2V_ 


+ AC 


“power “sting 


] * q * s * c 


Where: 


C - Baseline pitching moment coefficient as a function of 

“o o and M, (CMBASEH) . 

AC - Increment of pitching moment coefficient due to left and 

m 6_ right flap deflection as a function of o and M, (DCMFL 

f L,R and DCMFR) . 


AC 


m. 


- Increment of pitching moment coefficient due to left and 
right aileron deflection as a function of a, (HCHDAL and 


h,R HCMDAR) 


AC - Increment of pitching moment coefficient due to stabilator 

“STAB deflection as a function of a and M, (DCMST) . 

C - 3C /. , pitching moment coefficient due to pitch rate as 

m m dq 

q a function of M, CCMQH) . 



- 3 C / *, pitching moment coefficient due to angle of 
ni dot 

attack rate as a. function of M, CCMALPDIO . 


AC - Increment of pitching moment coefficient due to power 

“POWER effects as a function of a, Veq and 0^, (HCMPOW) . 

AC^ - Increment of pitching coefficient due to sting 

STING interference corra tion to wind tunnel data, (CMSTING) . 
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SF « 


[Cy * 6 + AC^ * 6 + (ACy + ACy ) + Cy * 6 R / 15 


R 


+ Cy * r * jr- ] q * S 
r T 


Where 

S 

A S 

AC., 


6 


- side force coefficient due to sideslip angle as 
a function of a and M, (HCYB) . 

- Increment of Cy due to flaps as a function of M, a 
and flap deflection, (HCYBF) . 

- Increment of side force coefficient due to left and 

l 6 right aileron deflection as a function of a, (HCYDAL 

A L,R and HCYDAR) . 

- 3Cy/8(5g, side force coefficient due to rudder deflection 

5 R as a function of M and a, (HCYDR) . 

- ^Cy/gj.* side force coefficient due to yaw rate as a 

r function of M, (HCYR) . 


o Rolling Moment 

RM-[C £ * B + AC £ * 6 + (AC £ + AC £ )+C *d R /15 

e 6 6 6 6- 

A L *r r 


+ (C* * r + C J *2vr + C 

r p T 


£ ] * q * S * b + Ac. 

POWER 6. 


Where: 

C. - 3C rolling moment coefficient due to sideslip 

l X. 

p angle as a function of M and a (HCLLB) . 




6 


- Increment of Co 


due to flaps as a function of M, 


and flap deflection, (RCLLBF) . 


a 


AC * 


6 




> 


R 


- Increment of rolling moment coefficient due to left 
and right aileron deflection as a function of a 
and M, (HCLDAR and HCLDAL) . 
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C 


£ s 


R 


• , rolling moment coefficient due to rudder 

deflection as a 'unctior. of c« and M, (HCLLDR) . 


C, 

x. 


T 


ac. 


AC, 


POWER 


• rolling moment coefficient due to yaw rate as 

a function of M and a, (HCLLR) . 

- SC^/.p, rolling moment coefficient due tc roll rate as 
a function of M and a, (HCMP) . 

- Increment of rolling moment due to asymmetric flaps , 
(HDCLLF) . 

- Increment of rolling moment coefficient due to power 
effects, CHCLLBF * BETA). 


o Yawing Moment 


YM- [C * 8 + (AC 
n 8 n 5 


+ dC 


) + ’ 


V 15 


+ (C * r + C 


p) * ] * q * s * b -i ic t 


Where: 


dC 




dC 


- 3C Q /jg, yawing momert coef f icier t due to sideslip angle 
as a function of M and a, (CNB) . 

- Increment of yawing moment coefficient due to aileron 
deflection as a function of M and a, (CNDAR and CNLAL) . 

- 8C /.. , yawing moment coefficient due to rudder deflection 

u o0«* 

as a function of M and a, (CNDR) . 

- 3C /. , yawing moment coefficient due to yaw rate as a 

XL oT 

function of M and a, (CNR) . 

* ?^ Q / gp . yawing moment coefficient due to roll rate as a 
function of a and M, (CNP) . 

- Increment of yawing moment due to asymmetric flops, 

(flCDNF) . 


3-22 



MDC A'910 
Volume 1 


QU 


Zxri 


3. 2. 3. 3 Center of Gravity Effects - The effects of center of gravity dis- 
placement are included in the aoment equations. All of the forces 

are referenced to a point on the aircraft at Fuselage Station 346.*, water- 
line c *. ?•: ion 96. 0 and coinciding with the X-axis of the aircraft. Incremental 
distances from the reference point to the center of gravity are computed and 
used in the following equations: 

o Pitching Moment 

PM - PM - F K * (346.6 - CGpsJ/ClZ) + F^ * (CG^ - 96.)/ (12) 
o Rolling Moment 

RM - RM - SF * (CG^ - 96.)/ (12) + F N * CG^/ '* 2) 
o hawing Moment 

YM = YM + SF * (CG p g - 346. 6)/ (12) + F A * C^/ (12) 

3. 2. 3. 4 Reference Axes Transformations - The aerodynamic force and moments 
are derived in the body axis system for use in the equations of notion. The 
following equations transform the high speed aero forces and mcmants fr >m the 
stability to body axis system. 


o L -Body Force (Positive down) 

F_ - -? T _ * COSq - ?__ A . * SINc 
Z_ LIFT DRAG 

o X-Bowv Force (Pocitive forward) 

% * - f drag * COSa + f lift * SINo 


o X-3ody Rolling Moment (Positive righ* - wing down) 

* Ki SIAB * C0S ° - ™S1U * SBi “ 

c Z-Body Yawing Moment (Positive nose right) 

* C0S ° + ^STAS * SINa 

3.2 4 ENGINE AND REACTION CONTROL SYSTEM FORjE AND MOMENT EQUATIONS - All 
*ngi^e and RCS forces and wonents are resolved in the oody a x is system 
(Figure 3.2-5) for use in *he equations of motion. The engine forces and 
moments include inlet momentum effects and scrub and boundary layer drag 
increment 0 . 


3.2. 4.1 Engine Threat Geomet ry - Figure 3.2-6 shows the relative geometry of 
the effective thrrat center with respect to the moment reference center used 
for the aerodynamic cU^ta. The following equations define the resulting engine 
thrust forces and moments about the \ 'rcraf t c.g. in the body axis system. 
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o X-Body Force (Positive forward) 

rr^ » t c * cos (©j + 1.5*) 

o Y-Body Force (Positive out right wing) 

rr_ » o. 

B 

o 2 ~Body Force (.Positive down) 
n z - -T c * SIN (!d J + 1.5*) 

B 

o Pitching Moment 


- T c * ((346.6 - TCFS)/12.) * SIS (0j + 1.5°) 

+ T c * ((96. - TCWL)/12.) * COS (Oj + 1.5°) - FT, * 
((CGFS - 346. 6) /12.) + FT- * ((CGWL - 96.)/12.) ® 


o Yawing Moment 

YHj - FT^ * CGBL/12 . + FXj * ((CGFS - 346.6) /12.) 
o Rolling Moment 

BMj - -FT * CBGL/12 . - FT^ * ((CGWL - 96.)/12.) 
o B 
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is the total thrust vector. 

TCFS is the thrust center fuselage station "ocation. 
TCWL is the thrust center waterline station location. 


P - 


Where: 


.. 

0F QUALITY 


3.2. 4.2 Ram Drag and Inlet Effects - Incremental forces and moments at the 
inlet are derived from ram drag forces as a function of inlet flow angle and 
from air mass flow as a function of aircraft bcay rates. The following 
equations summarize these inlet contributions at the aircraft c.g. The force 
application point on the inlet is defined as F.S, 220.63, W.L. 90.77 and B.L. 27. 


o X-Body Force 


r v 

o Y-Booy Force 


m 


Drag * COS 0 * COSa * q * AZ 

g 


m 


FI_ - -Ram Drag * Sing + -*p*AZ * r * AX 

g g 

o Z-Body Force 


FI_ ■ -Ram Drag * C0S£ * SINa + - * q * AX 

fl 8 

o Pitching Moment 

PMI * Ram Drag * COSB * ((STNa * AX) - (COSo * AZ)) 

- — * q * (AX 2 + AZ 2 ) 
g 

o Yawing Moment 

YMI « Ram Drag * ((COSB * COSa * CGBL/12.) - (SINS * AX)) 

+ — *p*AX*AZ- — *r* (AX 2 + AY 2 ) 
g g 

o Rolling Moment 

Rm - Ram Drag * ((SINS * AZ) - (COSB * SINa * CGBL/12.)) 

- — * p * (AZ 2 + AY 2 ) + -*r*AX*AY 

g g 


Where: 


Ram Drag * ^ * V^, lb^ 


AX - (f'GFS - 220.63) /12 . ' ft. 
AY « (27. - CGBL)/12. ~ ft. 
AZ - (CGWL - 90.77) /12. " ft. 
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3. 2. 4. 3 RCS Thrust Geometry - The relative geometry of the YAV-83 RCS is 
shown in Figure 3.2-7. To compute the pitch and roll forces due to the RCS 
requires accounting for the off-axis orientation of the pitch and roll valves. 
The following equations are used to derive the RCS forces and moments in the 
body axis system at the aircraft c.g.: 

o X-Bodv Force 

- FT rpv * SIN (7.63°) - FTjpy * SIN (8°) 

- « < 5 ‘> * < 7 *> * ("bo * FT KTOtT> 

NOTE: Subscripts defined in Figure 3.2-7. 

o Y-Bodv Force 

* "it - S » (! " > * 008 (7 *> * '"BW + rT STOD ) ET 
RCS 

(FT RVWD + ^RVWuW 



o Z-Body Force 

- FTppy * COS (7.63°) + FT Rpv * COS (8°) 
+ COS (5°) * COS (7«) *'(FT hwd + FT^) 


o Pitching Moment 

PMjgg » FT^ V * (COS (8°) * t 605 ---^'- 0 ^ + SIN (8°) * ( - - ) 

+ FT m * (COS (7.63°) * (^M; ggS) + si N (7.63°) * ) ) 

- COS (5®) * SIN (7®) * (" 3,8 j^ c S?fr) ) 

+ <"*Bn> + n m»D ) * ((C0S (5 ' ) * “* (7 *> * <- 32 ' ' S u: CGFS ) 

- COS (5®) * SIN (7®) * ( — ) ) 

o Y awing Moment 

YM RCS * ^TV * ( 12l ^ “ ((rr RVWD + ^RVWpRI 

- ^EVWD + *WlT> * < SIN (5 ^ * C0S (7 °> * 

+ COS (5°) * SIN (7 8 ) * ( 17 - 7 - 8 ^ - C - L ) ) 

+ (FT m * SIN (7.63°) - FT m * SIN (8°)) * CGBL/12 
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Yaw Valve (YV) 
FS 599.0 



(RVWU)rt 



AH dimensions in inches 
FS • Fuselage station 
WL - Water line 
BL • Butt line 


QP21-0M*4S 


YAV-88 REACTION CONTROL GEOMETRY 
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^•RCS * ((FT RVWD + ^RVWlpRT ~ (FT RVWD + ^NVWtpLT 5 * 

(SIN (5*) * COS (7°) * ( ^ L i2. - 83 ' * 89 ) + 008 (5 ^ * C0S (7 °' 

* C 17 ^; - (rr m * cos (7.63°) + rr m * COS (8°)) * 

CGBL/12 + FT^ * ( -° , - ^ CG ^ L ) 


3.2.5 WEIGHT AND BALANCE EQUATIONS “ The total weight of the aircraft is 
computed by summing aircraft dry weight (including the pilot), fuel weight, 
and water weight. Because the weight distribution in the aircraft is 
constantly changing due to fuel transfer and water usage, the center of 
weight moment for each factor is computed in all three body axes. In addition 
the moments of inertia and. the cross product of inertia, I_-, are computed as 
a function of these factors. 


The body axis coordinates of the center of gravity of each fuel tank and 
the water tank are stored in the subroutine, WTBAL07, data tables as a 
function of fuel and water remaining. The weight moment of the dry weight 
aircraft about the aircraft reference point in each axis is stored as veil. 
With the relative aircraft geometry determined, the body axis coordinates 
of the center of gravity are computed: 

- T«k (1) * ^ t ♦ I** «> * £ * ... 


\ 


WATER * 


2 0 


^2° 


X rr - (M_ + M_ + M_ )/WT , etc. 

*DRT yUEL *H 2 ° 


Similarly the y and z component of the center of gravity can be determined. 
The x component is referenced back to the wing leading edge and expressed as 
a percent of mean aerodynamic chord. 


CGPC - ((Xg. - 336.3)/ (99.792)) * 100 

The moments of inertia and product of inertia are computed in a similar 
manner using the parallel axis theorem: 

I n pJrl ■ T,I>k * (T Time 1 * **1 MU 1* + • • • 

^ ■ ttA:ER * 

■ Sry. + Sm 4 S 2 o 
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I X2 fuel " Tank (1) * ^TANK 1 * Z TANK l /r + * " 

hz* m * Am * V * W 8 
1x2 ■ Xn + Xll + ^20 

These Inertias are used along with the aircraft mass to determine the 
aircraft mass matrix which is used in the equations of motion. 

3.3 AIRCRAFT SUBROUTINES 


The YAV8B Program includes 6 subroutines which comprise the neces ry 
logic for math modeling the YAV-8B aircraft. The subroutines are the weights 
and balance routine (WTBAL07) , the primary flight controls routine (PFC07) , 
the reaction control system routine (RCS07) , the engine routine (ENG08) , 
the secondary flight control routine (SFC07) and the aerodynamics routine 
(AER0Y8B) . The sequential order of execution is shown in Figure 3.1-2. 


A brief description of the inputs r i outputs of each subroutine as 
well as top-level flow diagrams are pro' ded to familiarize the program user 
with the purpose of each subroutine. A more complete description of necessary 
program input and output variables is given in the user's guide, Section 3.4. 
The flow diagrams are simplified in order to aid the program user to interpret 
the program listings. Although the listings contain many useful comments, 
the flow diagrams can be useful if it \ acomes necessary to modify or change 
existing logic. «. more detailed description of the math model equations 
for the aerodynamic algorithms, the engine and RCS force and moment calcula- 
tions, and the weights and balance breakdown is given in Section 3.2. 


3.3.1 WTBAL07 SUBROUTINE - The WTBAL07 Subroutine (weights and balance) 
computes current values of the moments of inertia, center of gravity location, 
weight and mass of the YAV-8B. This subroutine interfaces with the engine 
subroutine (EG08) which determines the remaining internal fuel and water. 

The WTBAL07 subroutine is not designed to include stores configurations 
although it could be modified to simulate external stores cn the aircraft. 


The total weight of the aircraft is computed by summing aircraft dry 
weight, fuel weight and water weight. The moments of inertia about each 
body axis and the cross product of inertia, I^, are contf ruously computed 
as a function of the remaining fuel and water. The center of gravity 
location is determined by computing the weight moment of each factor (dry 
weifc-5 of aircraft, fuel weight and water weight) along the respective body 
axis and dividing by the gross weight. A selective bias in the weight 
moment for each axis is included if it is necessary to simulate a c.g. offset. 
This moment bias also will alter the respective moment of inertia. 


The inputs to this £’’ v t 
and lIDS positions. The 
weight and mass. The c f 
3.3-1 depicts the logi~ „'r 


fuel weight, water weight, and gear 
' vocation, inertias, gross 
*' the body axis system. Figure 
i • ■•/•^routine 


0 
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3*3.2 PFC07 SUBROUTINE - The YAV-8B primary flight control system consists 
of conventional ailerons, .^dder, and stabilator, with a Reaction Control 
System (RCS) which act c about all three axes during hover and transition. 
Longitudinal control of the YAV-8B is maintained by a combination of stabil- 
ator and reaction controls. Likewise, conventional ailerons, in combination 
with reaction controls, provide lateral control. Both longitudinal and 
lateral controls operate directly through mechanical linkage to the control 
stick. The directional control system is a combination of rudder and 
reaction controls as a function of rudder pedal position. The stabilator and 
ailerons can be ti immed by electro -mechanically moving the control stick. A 
limited authority Stability Augmentation System (SAS) provides rate stabil- 
ization. 

The PFC07 subroutine computes the control surface deflections and the 
RCS area openings. The control actuators are represented by first order 
transfer functions. The closed loop technique used for real time digital 
simulation of these actuators is based o^ equations derived from a pre- 
dictive method known as the "Reddy" metnod, (see Section 3. 2. 1.2). The 
math model includes rudder hinge moment limits applied to the total rudder 
travel , 

A simplified f ov diagram cf the PFC07 subroutine is shown in Figure 
3.3-2. Typical inputs include percent stick and rudder pedal, body rates, 
dynamic pressure and lateral acceleration. Outputs from the subroutine are 
surface deflections and RCS valve area openings. 

3.3.3 RCS 07 SUBROUTINE - The Reaction Control System (RCS) is used during 
hover and transition between full wingborne flight and fully jet -borne 
flight. The system is fully activated when the nozzles are rotated below 
15 degrees. It uses high pressure compressor bleed air, which is fed to the 
aircraft extremities and ejected through variable shutter valves connected 
to the aerodynamic surfaces (stabilator, rudder, etc.). 

The RCS07 subroutine computes the RCS pressure dynamics and mass flow 
characteristics. The RCS model in this subroutine is based on modified 
Rolls-Royce engine data and MCAIR determined aircraft installation losses. 
Included in the model is the RCS geometry for computing the resulting forces 
and moments in the body axis system (see Section 3.2.4). 

The inputs to the subroutine are the valve areas (calculated in the PFC07 
subroutine), compressor discharge temperature and pressure and nozzle deflec- 
tion. Other inputs include aircraft c.g. location and ambient pressure and 
temperature. The outputs are the total mass flow rate, the individual valve 
thrust values and the resulting forces and moments. The forces and moments 
are added to the engine produced forces and moments in the FNG08 subroutine. 
The forces and moments are in the body axis system. Figure 3.3-3 shows the 
logic x low used in the RCS07 subroutine. 
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FIGURE 3.3-2 

SUBROUTINE PFC07 FLOW DIAGRAM 


3-32 



' '■'T *1 ; - • MD A7S10 

Cr , 1 1 Volume I 

Ww “ v 



FIGURE 3.3-2 (CONTINUED) 
SUBROUTINE PFC07 FLOW DIAGRAM 
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FIGURE 3.3-2 (CONTINUED) 
SUBROUTINE PFC07 FLOW DIAGRAM 
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FIGURE 3,3-3 

SUBROUTINE RCS07 FLOW DIAGRAM 
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3.3.4 ENGQ8 SUBROUTINE - The ENG06 subroutine computes YAV-8B engine perfor- 
mance. The program data is based on Rolls-Royce engine data and mCAIR aircraft 
installation losses for an "average” YF4 02-RR-404 engine. Engine performance 
was obtained using a Rolls-Royce steady state engine performance computer 
program. Installation losses include inlet total pressure recovery, inlet 
drag, boundary layer bleed drag, compressor bleed (for aircraft services and 

the reaction control system) and mechanical power extraction. Engine performance 
is corrected for atmosphere c h a n ges from standard day. 

To give a representative engine, thrust, fuel flow, EGT, Pj, and fan RPM 
are assembled in data tables as functions of commanded throttle angle, air- 
craft Mach number, commanded nozzle angle and altitude. The steady state 
values are retrieved from the tables ar^ modified for the type of day, water 
and bleed corrections. The results are used to compute: 

o Engine forces and moments (including the RCS values computed in RCS07) 
o Reaction control system bleed 
o Fuel depletion 
o Water depletion 

Engine transients in RPM and EGT are modeled using first order time 
constants that vary with flight condition. The engine compressor discharge 
pressure limiter, EGT limiter and RPM limiter are also modeled. These are 
functional to limit fan speed and compressor discharge pressure. Selection 
of the limiters override causes all limiters except the compressor discharge 
pressure l i miter to be bypassed. 

There are a relatively large number of inputs tc the ENG08 subroutine. 

Some of the more Important inputs are percent throttle, altitude, atmosphere 
parameters, airspeed, c.g. location, initial fuel and water quantities, 
percent nozzle lever position and RCS forces and moments. Outputs include 
the total engine + RCS forces and moments In the body axis system, various 
engine performance parameters (RPM, EGT, etc.), nozzle position and re m a in ing 
fuel and water quantities. A simplified flow diagram of the ENG08 routine is 
shown in Figure 3.3-4. 

3.3.5 SFC07 SUBROUTINE - The YAV-8B secondary flight control system includes 

flaps and drooped ailerons. The flaps are operated by an electro-hyarsuiic 
system consisting of an electronic flap controller and two dual tandem 
actuators. Flap positioning is provided by the flap controller in accord^t- 
vith switch (mode) selection. There are two flap modes: up and down. The up 

mode provides either a constant 0 degree flap deflection or the capability to 
manually "beep” the flaps to any position from 0 to 25 degrees. The dovTi 
mode schedules flap deflection as a function of nozzle position and airste^. 

In addition the down mode switch also activates the aileron droop system, 
extending the ailerons to 15 degrees whenever the nozzles are 16 degrees or 
greater and the airspeed is less than 165 KIAS (with a 410 KIAS hysteresis) . 

Other secondary flight control systems modeled in the SFC07 subroutine 
include gear and LIDS deployment. Figure 3.3-5 outlines the SFCG7 subroutine. 

All extension and retraction rates , time delays surface load effects and 
essertial loeic characteristics relating to the YAV-8B secondary flight control 
systems a~e simulated. The inputs to the subroutine include nozzle position, 
flap mode selection, gear selection, LIDS selection switch, Leight-on-Wheels 
(WOW) switch and airspeed. Outputs are flap position, LIDS and gear positions 
and aileron droop position. 
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FIGURE 3.3-4 

SUBROUTINE ENG08 FLOW DIAGRAM 
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FIGURE 3.3-4 (CONTINUED) 
SUBROUTINE ENG08 FLOW OIAGRAM 
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FIGURE 3.3-4 (CONTINUED) 
SUBROUTINE ENG08 FLOW DIAGRAM 
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FIGURE 3.3-4 (CONTINUED) 
SUBROUTINE ENG08 FLOW DIAGRAM 
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FIGURE 3.3-5 

SUBROUTINE SFC07 FLOW DIAGRAM 
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FIGURE 3.:- 5 (CONTINUED) 
SUBROUTINE SFC07 FLOW DIAGRAM 
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FIGURE 3.3-5 (CONTINUED) 
SUBROUTINE SFCQ7 aOW DIAGRAM 
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3.3.6 AER0Y8B SUBROUTINE - The AER0Y9B Subroutine computes all of the aero- 
dynamic forces and moments acting on the simulated YAV-8B aircraft. These are 
computed through the use of stability and performance characteristics. The 
characteristics are stored in data tables and accessed through table look-up 
routines. When combined vith the aerodynamic reference dimensions of the 
YAV-8B (e.g., wing area, mean aerodynamic chord, and wing span) and the 
control surface positions, these coefficients are used to calculate total 
forces and moments in the body axis system. The aerodynamic equations for 
competing the total coefficients are given in Section 3.2.3. The aerodynamic 
data cover toe full range of aircraft performance, from hover to vingborne 
flight. The aerodynamic forces and moments are primarily functions of air- 
craft attitude, Mach number, power setting (with nozzles down) and control 
surface position. The baseline coefficients represent the aircraft with five 
bare pylons and the gun pods and LIDS onboard. Increments due to gear and 
LIDS extension, flap deflection and power effects are included. Ground 
effectc fcr altitudes from 0 to approximately 70 ft. are also included in the 
data base. The ground effects data include LIDS contributions. 

The AER0Y8B subroutine is divided into two distinct sections. The first 
section is the low speed (M <.0.3) region which includes ground effects. The 
aerodynamic data for this section are referenced in the body axis system. The 
second section is the high speed (M2 0*5) region. The data here are 
referenced in the stability axis system. For Mach numbers between 0.3 and 
0.3, table lookups are done in both sections and a linear interpolation is 
done to compute the respective total coefficient. The total aerodynamic 
forces and moments are deri d in the body axis system. 

The hover (u < 30 KTAS) and low speed regimes for the YAV-SB require 
special consideration when determining aerodynamic forces and moments (see 
Section 3. 2. 3.1). The low speed baseline aerodynamic coefficients are input 
for aircraft angles of attack between + 180 degrees. The angle of attack 
definition for the low speed table lookups is a -lend of body axis pitch 
attitude and true angle of attack. Likewise, in the lateral -directional 
axis the technique for determining aerodynamic forces and moments in hover 
differs from wingborae flight. For the aerodynamic table look-ups in hover 
(u < 30 KIAS) , the sideslip angle definition is a unique variable th't is a 
function of axial velocity, roll attitude and pitch attitude. For u > 50 RTAS, 
the conventional definition of 8 is used. For u from 30 to 50 RTAS a blend 
of the tr wo definitions is used. 

Power effects on the aerodynamic forces and moments are input as Incre- 
mental data. Oofc of the parameters for determining the power effects is the 
equivalent velocity ratio (Veq). This vari le is defined as ^li.O * Tf Thrust. 
This represent* the square root of the ratio of freestream dynamic pressure 
to the average nozzle exit dynamic pressure. 

Sene of the more significant inputs to the AER0Y8B subroutine are the 
control surface deflections, angle of attack, sideslip angle, gear and LIDS 
positions, body attitudes and rates, Mach number, velocity components in the 
body axis, dynamic pressure, thrust, nozzle position and flap position. 

Outputs include the total aerodynamic forces and moments in the body axis 
system and a large number of incremental aerodynamic coefficients. Figure 
3.3-6 represents the functional flow diagram for the AER0Y8B subroutine. 
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FIGURE 3.3-6 

SUBROUTINE AER0Y8B FLOW DIAGRAM 
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3.4 PROGRAM YAV8B USER'S GUIDE 


Execution of program YAV8B Is initiated by reading the EXECUTE deck ini o 
the computer. For most digital simulation runs, the EXECUTE, I50MS , and 
RTP T '."* card decks will contain the input variables required, depending on the 
c. . a. <•£ a desired run. Each of these decks is described in this section, 
are ti. 'able look-up subroutine, program defaults and overrides, array 
;i.rucvures . nd subroutine modification capability. 


3.4.1 VARIABLE ARRAY STRUCTURE - There are two primary arrays in the YAV8B 
prog x. for eouivalencing program variables. These two arrays are the F 
array (FARRAY) and the A array. The F array is dimensioned as F (3000) 
whereas the A array is dimensioned as A (1000). The first 1000 locations of 
the F array are generally reserved for non-aircraft dependent variables such 
as program directives. The first location of the A array is equivalenced to 
F (1001). Thus, A (100) is equivalent to F (1100), etc. This array is gener- 
ally used for aircraft dependent parameters such as aerodynamic coefficients, 
engine parameters, etc. 


Throughout the following discussion, the user will be referred to 
variables assigned to the F or A arrays. Program variables are equivalenced 
to F or A array locations for ease in assigning the variables common storage 
locations. Each subroutine of the program contains equivalence statements 
defining the variables required for input or output. A listing of the F and 
A array equivalenced variables is given in Appendix B of Volume II. The array 
locations that have no assignment are available as spares. 


3.4.2 EXECUTE DE CK - The EXECUTE deck loads files, reads date cards necessary 
for program initialization, executes program YAV8B, and executes the p-ogram 
SPRINT to output the desired parameters for a simulation run. Figure 3.4-1 
outlines the EXECUTE deck card sequence. 


3. 4. 2.1 Load and Execute Card Sequence The load and execute card sequence 
for the EXECUTE deck is shown in Figure 3.4-2. The executive first gets the 
needed compiled files, described below: 


FILE NAME 


DESCRIPTION 


YAV8B 

ISOMS 

YAER0 

YENGN 

YAVAC 

E0M 

ATMOS 

FMMPLY 

CARDS 

TBLKP 

RIP 

RTPDAIA 

MSGABT 

SPRINT 


MAIN PROGRAM; MUST BE LOADED FIRST 
TIMING SUBROUTINE 

CONTAINS SUBROUTINES AERODAT AND AER0Y8B 
CONTAINS SUBROUTINES YENGD AND ENG08 

CONTAINS SUBROUTINES AC07, RCS07 , PFC07, SFC07, WTBAL07 
EQUATIONS OF MOTION SUBROUTINES 

STANDARD ATMOSPHERE TABLE AND ATMOSPHERE SUBROUTINE ATM0S 
MATRIX MULTIPLY AND TRANSPOSE SUBROUTINES 
SUBROUTINE CARDS TO READ DATA CARDS IN EXECUTE DECK 

DATA TABLE LOOK-UP SUBROUTINES (FCALC, FSRCH) AND FUNCTIONS 
(FlA, FIB, etc.) 

SUBROUTINE RTF TO BUFFER DATA INTO STORAGE LOCATION 
SUBROUTINE RTPDATA TO SPECIFY WHICH DATA IS BUFFERED; 

MUST BE LOADED BEFORE RTP 
CONTAINS DIRECTIVES FOR PRINTING ERROR MESSAGES 
PROGRAM TO PRINT DATA TO OUTPUT 
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These are loader into temporary load file, YLEFL, after which the program 
is executed. During execution, the FTPDATA and RTP subroutines are buffering 
values of selected variables into temporary storage location TAPE3. Program 
SPRINT buffers these data out of storage and into a format suitable for the 
printer, based on a user's NAMELIST inputs (see section 3. 4. 2. 3). 

3. 4. 2. 2 Data Card Format - Program YAV8B requires initial values of some 
variables for execution. These data are entered through cards located behind 
the first end-of-record in the EXECUTE deck. The data cards are read by 
subroutine CARDS, which is called from program YAV8B. Card input format is 
shown in Figure 3.4-3. Column 5 contains one of the following control or 
data type identification numbers describing the data to be presented: 

COLUMN 5 INPUT INDICATION 

1 INTEGER DATA 

2 FI GATING POINT DATA 

3 OCTAL DATA 

4 ALPHANUMERIC DAxA 

5 CONTROL NUMBER INDICATING END OF A SET OF DATA CARDS 

6 CONTROL NUMBER INDICATING END OF ALL SETS OF DATA CARDS 

7 OPTIONAL ALPHANUMERIC DESCRIPTION OF TYPE 1 THROUGH 4 DATA; 

START 25 OR LESS CHARACTERS FOR COMMENTS IN COLUMNS 6, 

31, AND 56 

Up to three data points may be included on one data card for a given identi- 
lication number. Types 1 through 4 cards require array locations of the 
input variables (in the F array) , and each specification must end in columns 
10, 35, and 60. Input data itself is 20 or less characters ending in columns 
30, 55, and 80. 

Th-. inputs shown in Figure 3.4-3 are generally what is required for a 
simulation run. Only one data point per card has been input here for flexi- 
bility in changing cf single values. The data displayed is that used for a 
YAV-8B conventional Dutch Roll case. Variable names corresponding to the 
array locations of Figure 3.4-3 are defined below: 


F ARRAY LOCATION 

VARIABLE NAME 

DESCRIPTION 

150 

TST0PD 

LENGTH OF RUN IN SECONDS 

1011 

H 

ALTITUDE 

1030 

VT 

TRUE AIRSPEED IN FT/SEC 

1031 

GAMMA 

FLIGHT PATH ANGLE 

1040 

RMACH 

MACH NUMBER 

1103 

THETA 

PITCH ATTITUDE 

1117 

ALPHA 

ANGLE OF ATT-.CK 

1118 

BETA 

SIDESLIP ANGLE 

1195 

IN STD 

NON-STANDARD DAY FLAG 

1273 

BSM0MX 

INPUT FOR FUSELAGE STATION Cl 

1274 

BSM0MY 

INPUT FOR BUTT LINE CG BIAS 

1275 

BSM0MZ 

INPUT FOR WATER LINE CG BIAS 
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1495 

ISASL0N 

PITCH SAS SWITCH 

1496 

ISASLAT 

ROLL SAS SWITCH 

1497 

ISASDJR 

YAW SAS SWITCH 

1627 

IM0DE 

FLAP MODE SWITCH 

1635 

IFLAP 

MANUAL FLAP SWITCH 

1776 

PCTHR 

PERCENT THROTTLE 

1777 

PCN0Z 

PERCENT NOZZLE CONTROL 

1782 

PQFUEL 

FUEL WEIGHT 

1819 

PWATER 

WATER WEIGHT 

1882 

IGEAR 

GEAR UP /DOWN FLAG 

1959 

I CASE 

OPTIONAL CHECK CASE IDENTIF. ^ATION 


NUMBER 


3. 4. 2. 3 NAMELIST options - Program SPRINT is used to print data from the 
local file TAPE3. The printout consists of start of run data, time history 
data, and end of run data. The user may select various control options over 
the following categories via NAMELIST input: 

NAMELIST CONTROL CATEGORIES 

o Run Selection 
o Output Print Selection 
o Parameter Selection 
o Time to Print Selection 
o Names Override 
o Classification Override 
o Low Core Execution 

The data cards for the NAMELIST options are contained in the EXECUTE deck 
and follow the data cards described in the previous section. NAMELIST 
input options are available in each of these categories and are the following: 

Run Selection 

NRUNS Number of consecutive runs to bt printed, starting at 

the current file position. After NRUNS runs have been 
printed, the program will terminate wich a ST0P 1, 
unless a double end-of-file has been encountered before 
NRUNS has been satisfied. If this occurs, the program 
will terminate with a ST0P 2. Parameter NRUNS is not 
used i*. IFILEiKl) is non-zero. 

IFILES(20) File numbers of runs to be printed, starting at the 

current file position. File numbers must be in 
ascending order a no more than 20 separate runs may 
be specified by the user. A zero in the array termin- 
ates the sequence. After all runs have been printed, 
the program will terminate with a ST0P 1, unless a 
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IITLIS (20) double eno-of-file has been encountered before all runs 

(Cont'd) selected have been satisfied. ~f this occurs, the 

program Mil terminate with a STOP 2. Parameter IFILES 
overrides the parameter NEONS. 


Output Print Selection 


IBL0CK Format selection fn printing time history data. Legal 

values are as f .ol?c /s: 

Column Format (IBL3’2K“0) prints the complete time 
history of ten pan meters at a time in colums up to 441 
parameters out of *. maximum of 881 parameters . If the 
data exceeds 441 parameters it can be read by using 
NPLIM for the first 441 parameters and then using NPLIM 
in conjunction with the IXE array for the remaining 
parameters; see III array description. 

The Column T'crmat l squires more core since local files 
are used to s>rt the time history data for separate 
pages, and each file requires a separate buffer. If 
core is a problem the user should consider the Low Core 
Execution Optimum. 

Block Format (IBL(JCK“1) prints all time history para- 
meters for one time pass in a block, ten per line. The 
maximum number of parameters that can be printed is 881. 

Block Format wit a Names (IBL0CK»2) prints all time 
history paramete *s for one time pass in a block, five 
per line with pa amater names embedded within the block. 
The maximum number of parameters that can be printed is 
881. 

Column Format with Increased Accuracy (IBL0CK*3) prints 
the complete time history of five parameters at a time 
in columns up to 221 parameters out of a ma-srimmn of 881 
parameters. If the ciata exceeds 221 parameters the 
data can be read by using NPLIM for the first 221 para- 
meters and then using ITLIM in con j '-notion with the I7K 
array f or the remaining parameter-; see IXE array 
description for more details. 

This Column Format prints u. the values of the para- 
meters to 8 decimal places. And if used in conjunction 
with the IX0CT array "an print out the octal value of 
the parameter in 20 uutai digits (where 20 octal digits 
■ a. oG - hit word) . 
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Limit number of time hxbtcry parameters to be printed. 

If NPLIM is non-zero the first NPLIM parameters will be 
printed. Normally, all the time history’ pointers and 
names are printed on a start of run page. However, if 
parameter NPLIM is used, then only the pointers and 
names up to NPLIM will be printed. Also, for the 
Column Format printout, less core will be required, 
since only the sort file bufi >rs needed for NPLIM 
parameters will be allocated. 

Page print selection array. Each array element is a 
flag that corresponds to ter time history parameters 
for the Column and Block Format Options. For the 
Block Option with names each flag corresponds to five 
parameters up to a maximum of 881. For the Column 
Option with increased accuracy each flag corresponds to 
five parameters up to a maximum of 221. If IP flags 
are set to zero, the corresponding time history para- 
meters will not be printed, but the pointers and names 
will still appear on the start of run page. 

Parameter re-ordering array. For the Column Format 
(IBL0CK*O) this array allows printing up to 441 time 
history parameters in any desired order from a set of up 
to 881 time history parameters. For the Column Format 
(IBL0CK>3) , this array allows printing un to 221 time 
history parameters in any desired order from a set of 
up to 881 time history parameters. For both Block 
Formats this array allows printing up to 881 time 
history parameters in any desired order. Each array 
element corresponds to a parameter position. If any 
array element is non-zero, the original time history 
parameter corresponding to that element will be replaced 
by the parameter number indicated by the value of the 
array element. Only the parameters to be re-ordered 
need specification in the 1XR array. Also, for the 
Column Format Option some time history parameters can 
be duplicated to appear on more than one page, for 
example, by specifying the parameter number in more than 
one array position. This option can be used in combina- 
tion with the namelist parameters NPLIM and IP. 

NOTE : If the user wants Column Format (IBL0CK*O) and 

has a tape with NP *'442<NP<881) time history parameters, 
the first 441 parameters may be printed by using NPLIM* 
441. The remaining parameters (NP-441) may be printed 
by using NPLIM*NP-440 and IXR(l)*-442, which will auto- 
matically generate the correct IXR array. If the user 
wants Column Format with Increased Accuracy (IBL0CK*3) 
and has a tape with NP (222^<NP<881) time history para- 
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IXR(881) meters, a similar method is used to read the data in. 

(Cont’c) The firrt 221 parameters may be printed using N?LIM*221. 

If NP <_ 441, set NPLIM-NP-220 and IXR(l)— 222 and it 
will automatically generate the correct IXR array for 
this set of 221 parameters. 

If NP < 661. set NPLDf-NP-440 and IXR(l) —442. Other- 
wise, on this third pass leave NPLIM*22l, set IXR(l)— 422 
and it will generate the IXR array for parameters 442 to 
661. For printing out the remaining parameters from 
662<NP<881, set NPLIM-NP-550 and IXR(l) —662. 

IX0CT(881) Parameter octal print array. This array allows for 

printing of time history parameters in octal format. 

Each array element is a flag _hat corresponds to a time 
history parameter to be printed. That is, the second 
array element is associated with the second time history 
parameter to be printed, the third array element with 
the third parameter to be printed, and so on. This 
matchup is made regardless to whether or not the nth one 
to be printed is actually the nth t ime history para- 
meter or one that has been reordered using the IXR array. 
To have the parameter printed out in octal, its related 
XX0CT array element should be set high. 

For the Column Format with Increased Accuracy (IBL0CR-3) , 
which prints out five parameters per page, the local re- 
presentation is 20 characters long. For the Block Formats 
and the Column Format with ten parameters per page 
(IBL0CX»2,1,O) the octal representation is 10 characters 
long. The maximum msaber of parameters printed out for 
ISL0CK*O,1,2,3 is 441, 881, 881, and 221, respectively. 
These limits hold for the IX0CT array also. 

Time to Print Selection _ 


OP Print interval. Time history data will be printed once 

every data pass that TIME exceeds or equals (within a 
small tolerance) successive integral multiples of DP. 
Thus if the parameter DP is less than or equal to the 
smallest time increment, all the data will be printed. 

If DP*0.0, all passes of data will be printed regardless 
of the time increment This could be useful if , because 
of some error, TIME was equal to or less than a pre- 
vious pass. If DP»0.0, the namelist parameters TSTARI 
and TST0P are not operational. 
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TST0P (10) 


LAS TP AS 


Names Override 


INAMES 


Start /Stop times. Sections of time history data to be 
printed can be specified by use of the TSTART/TST0P 
arrays. The array elements, up to a limit of ten, 
correspond to the sections to be printed. The start 
time is specified by the TSTART element and the stop 
time is specified by the TST0P element. The stop time 
must be greater than or equal to the start time, and 
each start time must be greater than or equal to the 
previous stop time. All array elements beyond those 
desired must be zero. 

Print last pass of data. This parameter, when non-zero, 
causes the last pass of time history data to be printed 
regardless of the print interval or start /stop times. 


Names override. When this parameter is non-zero, the 
program will read data cards to input parameter names 
up to a maximum of 881. If names are on the data tape, 
they can be overridden using this option. If there are 
no names on the data tape, they can be supplied via data 
cards for the printout. Data cards must be in an (110, 
7A10) format immediately following the namelist cards. 
The first field is for an index indicating the para- 
meter number of the first name appearing in the second 
field. Up to seven names may be input per card. All 
the names on a card will be input in sequential order 
starting with the index for the first name. A blank in 
the name field will terminate input from that card. If 
seven names are put on a card, the names on the follow- 
ing card will be input sequentially even without index 
specification. A blank card terminates the data card 
input. 

If the program. is executed more than once in the same 
job, the names need to be entered only from the first 
namelist. If INAMES*2 , names entered from a previous 
nan ^ist will be de-selected, and new names can be 
^ ,ered by the user. 


Classification Override 


ICLASS Classification override. The classification is normally 

specified in the data tape start of run data, and will 
print out automatically unless this override option is 
specii^ed by the user, Legal values and the resulting 
printout are as follows: 

ICLASS*1 - Printout will have no classification, blanks 
will be printed instead of start of run classification 
data. 
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1CLaSS* 2 - Printout will have "Confidential" classifi- 
cation. 

ICLASS»3 - Printout will have "Secret" classification. 


Small buffer option for Column Formats printout. For 
these printouts (IBL0CK*O,3) , each sort file has a 
buffer which is defaulted to a length of 1002B for 
maximum efficiency. If the parameter L0WCORE is set 
non-zero, the buffer length is set to 102B, so there is 
a cors savings of 700B for every t'm parameters to be 
printed. The cost does not differ significantly, so 
the normal option is recommended because the job runs 
faster. The low core option can be used, however, when 
the available core is limited. This option does not 
apply to Block Formats (i.e. IBL0CK*1 or 2), which 
always use less core since no file buffers are necessary. 

Program SPRINT can be executed more than once in the same job in order 
to change the NAMELIST parameters 'for a different run or series of runs. All 
NAMELIST parameters not changed will retain their previously assigned values. 
The only exception is the parameter INAMES . This parameter will not stay 
non-zero on a subsequent execution so names do not have to be re-entered 
unless desired by the user. However, the names will be retained unless the 
parameter INAMES is set non-zero and new names are entered via data cards. 

Figure 3.4-4 shows the NAMELIST options input for the YAV8B conventional 
Dutch Roll simulation run. This set of cards is inserted behind the second 
end-of-record card in the EXECUTE deck. A summary of NAMELIST parameters 
is printed in the output for each run. and shown for this case in Figure 
3.4-5. The column format has been selected, and one page of the time history 
output for this run is shown in Figure 3.4-6. An example of block format 
output can be seen in Figure 3.4-7. 

3.4.3 I50MS DECK - Calls to the YAV8B aircraft subroutines and equations of 
motion are made in the I50MS subroutine. Within this subroutine, both 
dynamic time histories and static check cases can be modeled. In addition, 
initialization of variables can be performed and miscellaneous parameters 
can be defined. These features can be incorporated into the simulation with 
the use of the I50MS modify card deck. This deck modifies the ISOMS sub- 
routine to include any pilot control input time histories and/or program 
initializations . 

3. 4. 3.1 Dynamic Time Histories - Any simulated pilot inputs during a time 
history rur. are modeled in ISOMS before the operate loop call to AC07, a 
subroutine which calls the aircraft subroutines and equations of motion in 
the proper sequence. The pilot inputs can include stick and rudder pedal 
movements as well as throttle and nozzle changes. Any cockpit switches which 
need to be activated during a run, such as flap mode switch or water switch, 
can also be modeled. 


ICLASS 
(Cont ’d) 

Low Core Execution 
L0KCORE 
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The simuJation require, that stick inputs be made in percent of full 
stick deflection in both the longitudinal (100% full aft, 0% full forward) 
and lateral (-100% full left, 100% full right) axes. Rudder pedal is input 
in percent of total rudder pedal deflection (-100% full left pedal, 100% 
full right pedal). Varying throttle and nozzle control inputs are also made 
in percent for use in subroutine ENG08. See Figure 3.4-8 for commanded 
surface positions, R?M and nozzle angle as a function of percent stick, 
rudder, throttle and nozzle lever respectively. 

Pilot control inputs to the program are normally modeled as functions of 
run time. Figure 3.4-9 shows the location in 150MS where control inputs have 
been made for a YAV-8B wingbome Dutch Roll case. Data tables have been 
utilized in storing both the percent stick and rudder pedal movement as a 
function oi time. The format of the table look-up functions used to access 
these data tables is discussed in Section 3.4.5. It is acceptable to program 
step vtiputs of the controls as e function of time using IF statements in lieu 
uf data tables. For example, this is the procedure to use when modeling 
changes in cockpi * switch positions. However, the table look-up functions 
provide the capability of modeling more realistically the pilot control 
inputs due to linear interpolation between data points. 

3.4. 3. 2 Static Check Cases - When the output cf a simulation run is desired 
simply as a function of input and no control variations, the control input 
definitions can either be modeled as constants in the ISOMS subroutine or 
defined in the EXECUTE deck data cards. It may at some time be desirable 

to isolate a particular subroutine in order to examine data or logic. The 
analysis of a particular subroutine can be simplified by excluding the 
equations of motion and any aircraft subroutines irrelevant to the topic of 
interest. This modification can be made by replacing the call to AC07 with 
a call to the desired subroutine(s) in the I50MS subroutine. These sub- 
routines must be examined for necessary inputs, however, because s^oe 
variables used in the logic will no longer be provided values when the 
irrelevant subroutines are excluded. Additional parameters can be initialized 
through data card inputs to the EXECUTE deck. 

3.4. 3.3 Parameter Definition and Initialization Capabilities - The user may 
define additional parameters of interest in the ISOMS subroutine. Once 
equivalenced to available array locations, the parameters can be output in 
the RTPDATA subroutine (see Section 3.4.4). The advantage to defining a 
variable in ISOMS is that the definition can be made after calculations within 
all subroutines have been completed for a given pass. The user can be assured 
that the variable is using curre*.* and not past values in its definition. 

The entry point I5O0NCE is contained in 150MS , and provides an alternate 
means of initializing parameters for use in a simulation run. Greater flexi- 
bility is generally available in initializations made through the EXECUTE 
deck data cards. For necessary inputs which the user will not often be 
varying, however, entries made in I5O0NCE before the call to AC07 (' otne^ 
subroutines) can eliminate having to handle more EXECUTE data cards than are 
needed. 
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3.4.4 RTPDATA DECK - Program variables in common arrays can be deleted from 
or added to the output listing by modifying the subroutine RTPDATA with the 
modify deck RTPDATA. Subroutine RTPDATA specifies which array locations will 
be buffered out during a simulation run. The user needs to know the array 
locations of variables desired to be a part of the output. The format for a 
variable print-out request is: 

N»N+1 $ IB(N)-INa + xxx $ XN(N)-10H ccc, 

where "xxx" is the array location of the output, and "ccc" represents a ten 
character name which will identify the data of location " m " in the print- 
out listing. INA in this statement specifies that the data is coming from 
the A array, which begins at F(1001). If the variable address is defined in 
the P array and is less than 1000, INF should be used instead of INA. 

The capability exists for the user to output start of run data, time 
history data, and/or end of run data. Start of run and end of run variables 
are printed only once in the output, while time history variables can be 
recorded throughout the simulation run. The same statement format shown 
above is utilized for any of these three cases. Placement of a variable 
print-out statement within subroutine RTPDATA will determine the order that 
the variable name and its value are printed. An example of the output 
variables in RTPDATA specified for a YAV-8B wingborne Dutch Roll case is 
shown In Figure 3.4-10. In this example, time history and start of run variables 
will be displayed, but no end of run output will be produced. 

3.4.5 DATA TABLE LOOR-UP FORMAT - Data tables such as those containing aero- 
dynamic and engine performance data are accessed using subroutines and func- 
tions contained in file TBLKP . The appropriate interpolation ratio, index, 
and data point selections are based on given independent % . .■■'.able values and 
inputs provided in the TBLK? functions and subroutine calls. 

3. 4.5.1 Data Table Structure - The aerodynamic and engine data tables are 
defined in subroutines AER0DAI and YENGD respectively. The tables may be 
one-, two-, or three-dimensional. The order of the data within the tables 
is important when using a table look-up function. The first dimension of a 
data table is a function of the first independent variable, the second 
dimension likewise is a function of the second independent variable, etc. 

Figure 3.4-11 displays the structure for three-dimensional table CNFLAPT. 
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DIMENSION 2 INDEPENDENT VARIABLE - FLAP DEFLECTION 
DIMENSION 3 INDEPENDENT VARIABLE - NOZZLE ANGLE 


FIGURE 3.4,11 
DATA TABLE STRUCTURE 


3. 4. 5. 2 Fuse t Ions - The format for a one -dimensional table look-up with no 
previously calculated ratios or indices is: 

X0UT - F1A(RIN , NUMPTS , BOUNDS , VALUES 1 

where RIN is input independent variable 
NUMPTS is number of points 

BOUNDS is minimum ant n^ximum values of input variable, defined in 
subroutine DATA statements, and 

VALUE is table of outout function values corresponding to each breakpoint. 

For a two- c r three-dimensional table look-up, the arguments have the same 
meaning, but RIN will need to be specified as a temporary array whose elements 
are in two (or three), independent variables: 

TEMP(1)-RIN1 

TEMT(2)-RIN2 

X0UT-F2A ( TEMP , NUMPTS , BOUNDS , VALUES ) 
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Here also NUMPTS becomes a two-dimensional array defined in subroutine data 
statements an;* expressed as: DATA NUMPTS number in dimension 1, number in 

dimension 2 . BOUNDS becomes a four- rather than two-dimensional array, 
de^ini'tg the endpoints in both dimensions. 

Whenever a table look-up function is used, the interpolation ratios and 
indices for each dimension of the function are stored in a temporary location. 
If the next xogical function uses the same independent variable, number of 
points _id boundaries as arguments, the previously calculated ratios and 
indices can be retrieved from storage fbr use in this interpolation. Seven 
dif**rent functions are available in TBLKP: 


DIMENSION 

FUNCTION 

USE (RATIO MEANS RATIO AND INDEX) 

w*r (x) 

FlA 

Find x ratio; calculate f(x) 


FIB 

Calculate f(x) using previous x ratio 

w-f (x,y) 

F2A 

Find x,y ratios; calculate f(x»y) 


F2B 

Use previous x ratio; find new y ratio; 
calculate f(x,y) 


F2C 

Calculate f(x,v) using previous x,v ratios 

w*f (x,y ,z) 

F3A 

Find x,v ,2 ratios; calculate f(x,y,z) 


F3D 

Calculate f(x,v,z) using previous x,y,z 
ratios 


The breakpoints of all curves must be evenly spaced whenever an "A" 
function is used to calculate the ratios and indices. These functions do 
not extrapolate beyond the boundaries; chey interpolate within the boundaries 
only. Out-of-bounds input arguments will receive the closest boundary value 
of the function. Previous ratios can be calculated by an "A" or H B" function, 
or by subroutines FSRCH or FCALC. 

3. 4. 5. 3 Subroutine FSRCH - When the breakpoints of a curve are unevenly 
spaced, a call to FSRCH will initiate a calculation of the interpolation 
ratio and index to be stored and subsequently ’ised by one of the table look- 
up functions. Each dimension requires a unique call to FSRCH if more than 
one dimension has uneven breakpoints. The format for a call to FSRCH is: 

CALL FSRCH (RIN, BRKPT, NUMPTS , IPREV , IDIM) 

where RIN is input independent variable, 

BRKPT is array name of uneven breakpoints defined in DATA statement, 
NUMPTS is number of breakpoints in BRKPT array , 

IPREV is index position where search is to begin, a^d 
IDIM is for which dimension in the data table the* search is being done 
(1,2, or 3). 

The interpolation index f^r a given call to FSRCH is stored in IPREV, sc that 
this index will be a starting point the next time this particular FSRCH call 
is made. IPREA. must therefore be a unique name for all calls which do not 
have the same independent variable and breakpoint table specified in the call 
arguments . 
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3. 4. 5 . 4 Subroutine FCALC - FCALC assumes that breakpoints for a curve are 
evenly spaced. The format for a call to FCALC is: 

CALL FCALC (SIN, BOUNDS ,NOMPTS , IUIM) 

where SIN is input independent variable, 

BOUNDS is low and high boundary of input variable, 

NUMPTS is number of points in direction considered, and 
IDIM is for which dimension in the data tables the search is being 
done. 

An interpolation ratio and index are calculated and stored for use by one of 
the table look-up functions. 

3.4.6 PROGRAM DEFAULTS AND OVERSIDES - Each subroutine comprising the YAV8B 
program may have a number of variables assigned with default values . The 
program user should be familiar with the more important variables with 
default values. The following paragraphs describe these variables for the 
YAV8B, ATMOS and aircraft routines. Overrides to the defaults can be incor- 
porated in the data cards of the EXECUTE deck. 

3. 4. 6.1 YAV8B - All variables defined in the DATA statements can be changed, 
but the following will be of chief interest to the user: 

NPASSD Number of passes in reset; currently defaulted to 400 for 
dynamic stabilisation of the simulation. 

FDT Basic program update time increment in seconds: def atilt 

value is .05. 

CDT Current update time increment in seconds;, default value is 

.05. 

TSTOPD Run time limit; usually specified in the EXECUTE deck, but 
defaulted to 10 here in case user inadvertently excludes 
this input (in seconds) . 

DTBUF Rata at which subroutine RIP will buffer data to storage; 

should coincide with value of CDT (in seconds). 

3. 4. 6. 2 ATMOS - The basic ATMOS subroutine defines temperature, speed of 
sound, pressure, and density as a function of altitude for a standard day. 

If a non-standard day moael is desired, the user must modify the ATMOS stib- 
routine to include the values for these atmospheric variables In addition, 
a flag must be set to integer 1 in the EXECUTE data cards to bypass the 
standard day calculations in favor of the non-standard day inputs. This flag 
is INSTD, located in A(195) and specified as F(1195) on the data cards. 

3. 4. 6.3 AERPY8B - No formal options exist in the aerodynamic subroutine to 
affect the logical sequence -f computations. Any inputs necessary which are 
relative to cockpit controls cr aircraft configuration are provided by other 
subroutines. 
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3. 4. 6.4 ENGQ8 - The engine routine has a number of switches to consider, 
which are summarized below: 


F- ARRAY LOCATION VARIABLE NAME DEFAULT VALUE 


DESCRIPTION 


1781 

IDEFUEL 

1 

Flag for fuel bum during 
run; default value in- 
hibits fuel depletion 

1835 

IWATER 

0 

Water switch off 

1831 

IENGF0 

0 

Engine normal; no flame-out 

1874 

ILIM?FF 

« 

0 

All engine limiters on; 
when set to 1, EGT and TPM 
limiters are bypassed 

3. 4. 6. 5 RCS07 - No formal options exist in the reaction control system 
subroutine to affect the logical sequence of computations. Variable PTBLD 
is initialized to (1.) in I5Q0NCE because it is used in tnp denominator of an 
equation here before it is defined by subroutine ENG08. 

3.4. 6.6 PFC07 - 
a value of zero. 

Stability augmentation system switches are defaulted off with 
The variables and locations are: 


PITCH SAS 
ROLL SAS 
YAW SAS 

ISASL0N 

1SASLAT 

1SASDIR 

F(1495) 
P (1496) 
F(1497) 

Although it is not necessary to include these flags as inputs in the EXECUTE 
data cards when the switches are off, they are included simply for user 
information. 

3. 4. 6. 7 SFC07 - The secondary flight controls subroutine has a number of 
options available which are summarized below: 

F-ARRAY LOCATION 

VARIABLE NAME 

VALUE 

DESCRIPTION 

1882 

J.GEAR 

0 

1 

Default; gear up 
Gear down 

1628 

ILID 

0 

1 

Default; LIDS will extend and 
retract with gear 
Emergency LIDS retract 

1627 

IM0DE 

0 

Default; flaps in UP mode; move- 


ment of flaps must be commanded 
by IFLAP switch 


1 Down mode; flaps operate automat- 
ically as a function of nozzle 
position and airspeed; aileron 
droop activated in this mode; 
IFLAP ignored 
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F- ARRAY LOCATION 

VARIABLE NAME 

VALUE 

DESCRIPTION 

1635 

IFLAP 

0 

Default; no manual flap movement 



-1 

Retract flaps as long as IFLAP 
has this value 



1 

Extend flaps, as long as IFLAP 
has this value, up to 25° 

1900 

IW0NW 

0 

Default; no weight on wheels 



1 

Weight on wheels; should be set 
when aircraft is on the ground 


If the flap mode switch is J? and the user wishes to extend the flaps less 
than 25°, logic must be inserted into ISOMS, before the call to AC07, to 
limit the duration of time that the IFLAP switch is set to 1. 

3. 4. 6. 7 WTBALQ7 - A user may wish to bias the center of gravity to match a 
particular test configuration. This can be done by adding increments to the 
tctnl moments used in the cer>te:r of gravity calculations. The defaults for 
all biases are zero. 

When a particular C.G. is desired, the total moment needed can be cal- 
culated for a given aircraft weight, i.e. , for the fuselage station C.G.: 

T0TM0MX«CGFS*WE1GHT 

where T0TM0MX-DRYM0MX + FUELMX + . . . + BSM0MX 

Since the moment contributions due to the dry aircraft (DRYM0MX) , fuel 
(FUELMX) , etc., and the total aircraft weight are known , the bias (BSM0MX) 
required for a particular C.G. locat .on may be determined. The same procedure 
is followed for butt line and waterline C.G. biases, and the results for each 
are entered on EXECUTE data cards in the following locations: 

BSM0MX FC1273) 

BSM0MY F(1274) 

BXM0MZ FC1275) 

These biases are added into the total moment calculations used to calculate 
the aircraft C.G. location in WTBAL07. 

3.4.7 SUBROUTINE MODIFICATION - M0DIFY is a utility available on the Network 
Operating System (NOS) for Control Data Corporation (CDC) Cyber 170 Series 
Computer Systems; CDC Cyber 70 Series, Models 71, 72, 73, and 74 Computer 
Systems; and the CDC 6000 Series Computer Systems. M0DIFY is used to main- 
tain and update subroutines that are on the decks of modify file YAVM0D. 

These decks and the subroutines defined in them are summarized below: 
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QUALITY 


DECK OK YAVM0D SUBROUTINES IK DECK 

YAV8B (Program) YAV8B 

I50MS ISOMS 

YAER0 AERO DAT 

AER0Y8B 

YENGN YENGD 

ENG08 

YAVAC AC07 

RCS07 

PFC07 

SFC07 

WTBAL07 

ATMOS STNDAY 

ATMOS 

RTPDATA RTPDATA 

Figure 3.4-12 is an example of a modify card deck for incorporating 
changes into subroutine ISOMS, Insert and delete cards are placed between 
the *DECK and *EDIT cards after the first end-of-record . The insert card 
will be of the format (beginning in column 1) : 

^INSERT or *1 nrrn 


where nnn is a deck line number. Deck line numbers are those displayed to 
the right of subroutine statements, as shown in the listings of Volume II, 
Appendix A. Input cards following a given insert card will be placed by 
MODIFY in the ISOMS deck immediately after nnn. The delete card will have 
the format (beginning in column 1) : 

* DELETE or *D nnn 

or *D nnn, mmm 

where nnn, mm m specifies the first ani last deck line numbers of a group of 
lines to be deleted. Additional statements on cards following a given delete 
card are permitted, and these will replace the lines removed. 

Tc make modifications to a deck other than I51KS in Figure 3.4-12, the 
user must replace the name I30MS on the FTN, REPLACE, *DECK, and *EDIT 
cards with that of the desired deck. Further details of the M0DIFY utility 
can be found in Reference (7). 
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FIGURE 3.4-12 

MODIFY CARD DECK FOR I50MS 
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4. YAV-8B LINEAR MATHEMATICAL MODEL 


This section presents the linear mathematical model of the YAV-8B air- 
craft. The model is composed of small perturbation stability derivatives in 
the form traditionally used in the analysis of conventional aircraft. 

AV-8 project experience ha^ shown that this form is adequate for control 
system design and handling qualities analysis in the V/STOL flight regime 
providing the linear analysis is supplemented by nonlinear six degree of 
freedom analyses. The nonlinear analyses can be obtained through real time 
pilot in the loop flight simulation or non-real time digital computer runs. 

4. 1 LINEARIZED EQUATIONS OF MOTION 

The linearized equations of motion in stability axes are presented in 
Figure 4.1-1. The corresponding sign convention and nomenclature are pre- 
sented in Figures 4.1-2 and 4.1-3. The variables in the equations are 
perturbation quantities from the initial condition. 

4.2 TRIM CONDITIONS AND STABILITY DERIVATIVES 


The stability derivatives presented in this section include the effects 
of aerodynamics, gross thrust and inlet momentum. These derivatives were 
obtained by perturbing the nonlinear model about a 3 degree of freedom wings 
level trim condition. Using a 3 degree of freedom trim ensures that the 
effectiveness of the RCS, which is a strong function of engine fan speed, is 
represented accurately in forming the control effectiveness derivatives. 
Figures 4.2-1 and 4.2-2 present the trim conditions and the corresponding 
stability derivatives used to generate the linear results that are compared 
to YAV-8B flight test data in Section 5. The corresponding weight and 
balance data are presented in Figure 5.1-5. Figures 4.2-3 to 4.2-5 present 
trim conditions and the corresponding stability derivatives for 24 flight 
conditions that characterize the YAV-8B throughout its flight envelope. The 
weight and balance data are presented in Section 2.3. 
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VARIABLE 

DEFINITION 

UNITS 

A NT, 

U, L. R 

Aircraft nose down, up, left, right 

• 

Z 


Mean aerodynamic chord 

in 

FRL 


Fusalaga Reference Lina 

- 

FS 


Fuselage Station 

in 

I 


acceleration due to gravity 

32.174 ft/ 

G*ar 


Landing Gear 

- 

h 


Altitude 

ft 

I 

XX 

iyy 


Moment of Inertia about X-axis 
Moment of Inertia about Y-axis 

Slug-ft 2 

Slug-ft 2 

I zz 


Moment of Inertia about Z-axis 

Slug-ft 2 

I TZ 


XZ Product of Inertia 

Slug-ft 2 

L 


Rolling Moment (positive is KWD) or Left 

ft-lb, - 

LID 


Lift Improvement Device 

- 

M 


Pitching Moment (positive is ANU) 

ft-lb 

N 


Yawing Moment (positive is ANR) 

ft-lb 

n f 


Engine Fan Speed 

l 

n 

z 


Load Factor 

- 

P 


Rolling Velocity 

rad/ sec 

q 


Pitching Velocity 

rad/ sec 

q 


Dynamic Pressure 

psf 

r 


Yawing Velocity 

rad /sec 

RCS 


Reaction Control System 

- 

RWD 


Right wing down 

- 

TEU t 

D, L. R 

Trailing Edge Up, Down, Left, Right 

- 

u 


Forwarc speed 

* '/see 

V 


Total freestream velocity 

ft/sec 

tfL 


Waterline 

in 



Force in X direction (positive is forward) 

- 

Y 


Fore* in Y direction (po*itiv« i* out right wing) 

Z 


Force in Z direction (positive is down) 

- 

a 


Fuselage (WL) angle of attack 

radian 

6 


Angle of sideslip 

radian 

r 

0 


Initial flight path angle 

radian 


FIGURE 4.1-3 

LINEARIZED MATHEMATICAL MODEL AND STABILITY DERIVATIVES 

NOMENCLATURE 
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VARIABLE 

DEFINITION 

UNITS 

6 a 

Aileron deflection (positive is left aileron 
TED' 1 

radian 

6 ad 

Drooped aileron deflection (positive is TED) 

radian 

6 arcs 

Aft RCS deflection in equivalent horizontal 
tail deflection for 6g >_ 2° ^ARCS * 

6 R < 2 ° f ARCS “ 0 

radian 

F 

Trailing edge flap deflection (positive is TED) 

deg 

6 frcs 

Forward RCS deflection in equivalent horizontal 
tail deflection for 6 jj >_ 2° SjrcS * ® 

< 2 ° 5 FRCS * 6 H 

radian 

°H 

Horizontal tail deflection (positive is TET) 

radian 

5 r 

Rudder deflection (positive is TEL) 

radian 

6 rrcs 

Roll RCS deflection in equivalent aileron 
deflection & %RCS “ °A LEFT 

radian 

6 yrcs 

Yaw RCS deflection in equivalent rudder 
deflection ^YRCS " S R 

radian 


Roll angle (positive right w ^g down) 

radian 

e 

Pitch angle 

radian 

V 

Yaw angle 

radian 


FIGURE 4.1-3 

LINEARIZED MATHEMATICAL MODEL AND STABILITY DERIVATIVES 
NOMENCLATURE (CONTINUED) 
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M * 

u 

CoSVc/Iyy)^ 

*■ u 

(l/ffsec) 

L s- 

(pSV^IaaOC. 

l a 

(l/sec 2 ) 

*h* 

(pSVH/2Iy7)C^ 

(l/ffs«c 2 ) 

l b’ 

(pSVb 2 /4Ixx)C. 

6 

(l/sec) 

M * 
a 

(pSV^/ZIyyK 

ul 

a 

(1/sec 2 ) 

L r " 

(pSVb 2 /4lxx)C 1 

r 

(l/sec) 

m 

M* 

a 

(pSVc 2 /4Iyy)C 

IB* 

a 

(1/sec) 

V 

(pSVb 2 /4Ixx)C t 

P 

(l/sec) 

M * 

q 

(pS7c 2 /4l77)C m 

(1/sec) 

L 6 * 

(pSV 2 !,/ 2Irt)C # 

(l/sec 2 ) 


(pSV^c/ZIyyJC 

m 6 

(1/sec 2 ) 

V 

(pSV^/ZIZztt 

*8 

(l/sec 2 ) 

X « 

U 

(pSV/m) C x 
u 

(1/sec) 

n b* 

(pSVb 2 /4Izz)C 

n 6 

(l/sec) 


(pSV/2m) C_ 

(1/sec 2 ) 

N - 

r 

(pSVb 2 /4Izz)C 

°r 

(l/sec) 

x - 

a 

(pSV 2 /2m) Cj 
a 

(f t/sec 2 ) 

V 

(pSVb 2 /4Izz)C 

s 

P 

(l/sec) 

X* • 
a 

(pSVc /4a) 

0 

(ft/ sec) 

N 5- 

(pSV 2 b/2Iz2)C 

n 6 

(l/sec 2 ) 

X - 

q 

(pSVc/4m) C_ 

q 

(ft/sec) 

T s- 

(pSV i /2m) C 

y B 

(f t/sec 2 ) 


(pSV 2 /2m) C_ 

h 

(ft/sec 2 ) 

Y r 

(pSVb/4m) C 

7 8 

(ft/sec) 

2 ■ 

u 

(pSV/o) C 2 

U 

(l/sec) 

T • 

r 

(pSVb/4m) C 

7 r 

(f t/sec) 

2 h- 

(pSV/2m) C 

Ti 

(1/sec 2 ) 

Y - 

P 

(pSV b/4o) C 

7 P 

(f t/sec) 

2 a- 

(pSV 2 /2a) C 2 

a 

(ft/sec 2 ) 

Y 4* 

(pSV 2 /2m) C 

y 6 

(ft/sec 2 ) 

2* “ 
a 

(pSVc/4m) C 2> 
a 

(ft/sec) 




2 - 

q 

(pSVc/4m) C 2 

q 

(f t/sec) 




2 *‘ 

(pSV 2 /2m) C 

(ft/sec 2 ) 





FIGURE 4.1-3 

IINEARIZEO MATHEMATICAL MODEL AND STABILITY DERIVATIVES 
NOMENCLATURE (CONTINUED) 
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FIGURE 4.2-1 

TRIM CONDITIONS FOR TIME HISTORIES OF SECTION 5 
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LATERAL -DIRECTIONAL STABILITY DERIVATIVES FOR 
TIME HISTORIES OF SECTION 5 
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FIGURE 4.2-4 
TRIM CONDITIONS 
FULL FLIGHT ENVELOPE 
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FIGURE 4 . 2-4 (CONTINUED) 
TRIM CONDITIONS 
FULL FLIGHT ENVELOPE 
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FIGURE 4.:-4 (CONTINUED) 
TRIM COITIONS 
FULL FLJrj ENVELOPE 
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FIGURE 4.2-5 

LONGITUDINAL STABILITY DERIVATIVES 
FULL FLIGHT ENVELOPE 
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FIGURE 4.2-5 (CONTINUED) 
LATERAL-DIRECTIONAL STABILITY DERIVATIVES 
FULL FLIGHT ENVELOPE 
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FIGURE 4.2-5 (CONTINUED) 
LONGITUDINAL STABILITY DERIVATIVES 
FULL FLIGHT ENVELOPE 
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FIGURE 4.2-5 (PCNTINUED) 
LATERAL-DIRECTIONAL STABILITY DERIVATIVES 
FULL FLIGHT ENVELOPE 
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FIGURE 4.2-5 (CONTINUED) 
LATERAL-DIRECTIONAL STABILITY DERIVATIVES 
FULL FLIGHT ENVELOPE 
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5. i BACKGROUND 

The YAV-8B mathematical model is based on AV-8A and YAV-8B wind tunnel 
test and flight test results. The model has been periodically refined since 
the first AV-8A air-to-ground attack simulations at MCAIR in 1972. The 
development of the V/STOL model commenced in 1974 with the definition of an 
engineering simulation of a hover, short takeoff and landing and transition 
mode for the AV-8A. During this period the coordinate system used for hover 
and £he jet induced velocity parameter (Veq) were developed. Subsequent studies 

added data to simulate ground effects. To help validate the AV-8A simulation, 
John Farley, Chief Test Pilot for Hawker-Siddley , flew the simulator in hover 
using the motion base simulator. During 1975, additional improvements were 
made to the model including new air-to-ground weapon delivery modes and dis- 
plays, improved ship and land-based landing aids, and continued use of the 
model for air-to-air combat. All the changes were evaluated by the USMC 
using the flight simulation facilities. Between early 1972 and late 1975, 
there were 13 contractual AV-8A simulations. AV-8A modeling was refined 
throughout each of these programs based upon additional test data and 
engineering analyses. This simulation capability and experience permitted 
MCAIR to design and conduct a preliminary AV-8B Simulation Evaluation in 
June, 1S75. 

This evaluation utilized a previously-developed simulation of the AV-8A 
aircraft modified with AV-8B lift, drag and pitching moment characteristics. 
Following this evaluation, the AV-8B simulation was further modified to 
include wind tunnel-derived AV-8B aerodynamic and flight control characteris- 
tics in conjunction with AV-8B weight and balance, and propulsion character- 
istics. During October and December 1975, a second simulation was conducted 
to demonstrate flying qualities and perrcrmance improvements estimated for 
the AV-8B. Two experienced USMC AV-8A pilots participated in this simulation 
and their comments and evaluation contributed to improvements in the refined 
AV-8B model. 

Simulator modeling of the prototype AV-8B (YAV-8B) began in 1977. 

During 1978, the simulation was used for V/STOL flying quality evaluations 
and for pilot familiarization of the flight envelope prior to first flight 
of the YAV-8B. 

The YA T ' ^6 simulation was used very effectively during the entire 
prototyp ight test program in the following areas of activity: 

o Pilot familiarization prior the first flights 

o Determination of control time histories and flight profiles 
for maximum performance VTO's and STO's 

o Determination of handling qualities and pitch transients 
during landing approaches and transitions to and from 
wingbome and jetbome flight 
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o Simulated control failures 

o Flight trajectories following simulated engine flame-outs 

o Handling qualities in vingborne flight 

The YAV-8B simulator proved to be representative of the aircraft and very 
useful duri n g each of these activities. 

Prior to the YAV-8B engine/inlet compatibility testing in 1980, flight 
trajectories and approach/ landing techniques following a simulated engine 
flame-out (SFO) were developed and practiced on the simulator. Landing pro- 
files, flap angle, and the no 2 zle/pover setting needed to represent engine 
vindmilling were developed. During actual flight SFO practice, pilots 
reported that the flight SF0 ! s were very close to those experienced in the 
YAV-8B simulator. Also, the aircraft settings (idle pover/40° nozzle) chosen 
to simulate engine out conditions were representative of the actual engine 
out characteristics. 

The YAV-8B mathematical model is based on extensive flight test and 
wind tunnel test results. Figure 5.1-1 summarizes the flight test conditions 
flown. A total of 364 flights were flown prior to the start of a store 
separation flight test program in 1982. A summary of both aerodynamic and 
propulsion wind tunnel test programs is presented in Figure 5,1-2. The 
aerodynamic tests consisted of both powered and unpowered tests on a 15% 
scale and full scale models over the entire Mach range. They covered the 
jetbome, semi-j etbome and wingborne inodes of flight in and out of ground 
effect. Over 6600 hours cf tests are shown. A similar wind tunnel test 
program was performed on the YAV-8B inlet to provide a firm data base for 
the propulsion performance package. Approximately 2000 hours of wind tunnel 
testing are shown. 

The YAV-8B aircraft is equipped with both a noseboom angle of attack 
(AQA) probe and a production angle of attack indicator. The relationship 
between the noseboom ADA and the production AO A, in units, is shown in 
Figure 5.1-3. Note that the data fairing 

a FRL © * T * - 864 Production +1 

was used exclusively iu this report to compare flight test and predicted 
angle of attack. Previously published reports such as Reference 1 used 
production angle of attack exclusively. A similar situation exists with 
respect to the pitot static system position error corrections shown in 
Figure 5,1-4. Reference 1 used the estimated values while this report 
uses a fairing approximately halfway between the estimated and tower 
flyby results to compute calibrated airspeed. 
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FIGURE 5.1-1 

YAV-8B FLIGHT TEST ENVELOPE 
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FIGURE 5.1-2 

YAV-8B RELATED WIND TUNNEL TESTS 














A Mach Correction A Vc Correction KIAS A Hpc Correction ft 


Qf f 0^*' ^ 


Tower Flyby Data 
YAV-8B No. 1 Flight 42 

Symbol Gear Nozzle Fiao 


O 

□ 


Up 

Down 


Aft 

VAR 


5* 

VAR 


MDC A7910 
Volume I 





Note: A Corrections s True - Indicated dri4w« 

FIGURE 5.1*4 

YAV-83 PILOT STATIC SYSTEM ESTIMATED POSITION ERROR CORRECTIONS 
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The time history comparisons, presented in this section between flight 
test results and predictions from the YAV-8B mathematical models, were 
done using the weight and balance data presented in Section 2*3 adjusted to 
the actual weigh*: and c.g. location defined by the flight data. These weight 
and balance data are presented in Figure 5.1-5. 

5.2 JET BORNE AND SEMI-JET BORNE VALIDATION 


A time history comparison of predicted and flight test results during 
hover stability tests is presented in Figures 5.2-1 to 5.2-3. The predicted 
include both nonlinear and linear mathematical models defined previously. 

The predicted landing approach characteristics of Figure 5.2-4 shew 
good correlation with flight test results. The predicted data were obtained 
by fixing ^czzle angle and angle of attack and iterating on tail deflection, 
fan speed and airspeed until a 3 degree of freedom wings level trim was 
obtained. 

A comparison of angle of attack stability for flap settings of 25° and 
61.7° (STO) is presented in Figures 5*2-5 and 5.2-6. The predicted longitudinal 
stability and trim tail deflections are representative of the actual aircraft. 

Longitudinal dynamics in Semi-Jet Borne flight are compared in Figure 

5.2- 7 for both the nonlinear and linear mathematical models. 

The predicted Semi-Jet Borne SAS on steady heading sideslip character- 
istics shown in Figure 5.2-8 were obtained by using stability derivatives 
and a three degree of freedom approximation at the flight test condition. 

The predicted Semi-Jet Borne SAS on Dutch Roll characteristics are 
compared to flight test data and the requirements of MIL-F-83300 in Figure 

5.2- 9. The Level 1 requirements are satisfied. A time history comparison 
of Dutch Roll characteristics is presented in ~igure 5.2-10. 

5 . 3 WINGBQRNE VALIDATION 


A comparison of horizontal tail and angje of attack to trim at medium 
and high altitudes is presanted in Figure 5.T-1. 1 lie predicted data were 

computed using the same Mach number, nozzle angle and fan speed as the flight 
test data. 

The predicted and flight test short period response characteristics are 
compared in Figure 5.3-2 to 5.3-3. The short peric * frequency of Figure 

5.3-2 shows the characteristics are within the Level 1 boundary. The damping 
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FIGURE 5.2-3 

YAV-8B RESPONSE TO PEDAL INPUT IN HOVER 
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FIGURE 5. '.-4 

YAV-8B LANDING APPROCH CHARACTERISTICS 
STOL Flap, Ailerons Drooped 15* 
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Flight Data 

Sym 6j (deg) C.G. (%c) N f /V 0(%) V KCAS 6 F (deg) . ‘titude (ft) T(*F) 

□ 52 12.2 79 100 25 1.600 64 

A 63 13.2 66 130 25 3.000 59 



FIGURE 5.2-5 

YAV-8B SEMI-JET BORNE ANGLE-OF V. TACK STABILITY 

25° Flap Ailerons Not Drooped 
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YAV-8B SEMI-JET BORNE ANGLEOF-ATTACK STABILITY 
STO Flap Ailerons Drooped 15* 
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FIGURE 5.2-fi 

YAV-8B V/STOL STEADY HEAOING SIDESLIP CHARACTERISTICS 
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FIGURE 5.2*9 

YAV-8B SAS ON DUTCH ROLL CHARACTERISTICS IN SEMIJETBORNE FLIGHT 
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FIGURE 5.2-10 

YAV-8B RESPONSE TO RUDDER IN SEMIJET BORNE FLIGHT 

Gross Weight = 17,070 lb Nip = 88% 0j=61* 6 + = 25 ' 

Mach 0.24 Altitude = 4,700 ft 
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Flight Data _ 

Sym GW(lb) Alt(tt) N f l\6 [%) 

O 1 7700-1 8.78C 17,650-19.810 70.2-91.8 

□ 17,270-19,580 34,930-40,730 90.2-102.8 
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FIGURE 5.3-1 

YAV-6B WINGBORNE TRIM CHARACTERISTICS 
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FIGURE 5.3-2 

YAV-8B SHORT PERIOD FREQUENCY REQUIREMENTS 
Wingbome Flight 
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ratios of Figure 5 .3-3 are essentially Level 1 below 20,000 feet 
and Level 2 above 20,000 feet. However, MIL-F-8785B states that above 20,000 
feet the Level 1 damping requirements may be relaxed to the Level 2 require- 
ments. The longitudinal maneuvering characteristics are compared for the 
symmetrical pullup showr in Figure 5.3-4. The predicted characteristics 
correlate well with the flight data. 

Figure 5.3-5 summarizes roll rate capability and maximum sideslip angles 
generated during 1% rolls, ine flight data are primarily for low to moderate 
altitudes. The predicted characteris tics were computed for 7500 feet altitude 
using a nominal gross weight, and nominal command input time history. 

The comparison of steady heading sideslip characteristics presented in 
Figure 5.3-6 indicate the predicted static lateraJ^ directional stability 
and control characteristics are very representative of the actual aircraft. 

The predicted Dutch Roll charac teris tics are compared to flight test data 
and the requirements of MIL-F-8785B in Figure 5,3-7. The YAV-8B meets the 
Level 1 requirements below 20,000 feet and Level 2 above 20,000 feet. A 
time his ton’ comparison of Dutch Roll Characteris tics is presented in 
Figure 5.3-8. The predicted characteristics compare well with flight test 
resiT . 

5.4 PROPULSION SYSTEM VALIDATION 


The validity of the propulsion system model is demonstrated by the 
comparison of flight data and predicted performance in Figures 5.4^1 and 
5.4-2. These flight data include wet and dry Vertical Takeoffs (VTO’s) with 
stick neutral (minimum RCS reaction control system bleed), hence these cases 
are for a singular demand condition. The predicted engine acceleration 
characteristics show similar slopes as that of the flight data. Absolute values 
will not necessarily agree due to the following: 

o Engine to Engine Variation: The propulsion system model represents 

an ’'average” YF402-ER-404 engine. Actual engines used in the YAV-8B 
exhibited higher Exhaust Gas Temperatures and higher fuel flow at 
steady state conditions than an "average" engine. 

o Reaction Control Bleed: Bleed, which affects all engine parameters 

shown, is a function of control requirements which were not 
included in these cases. 

o Fuel Control Unit (FCU) Variations: Engine acceleration character- 

istics are influenced by FCU settings which differ with each 
engine/FCU comb ination. 

The stepwise characteristics of the measured fuel flow is due to the 
low sampling rate of the flight data system and does not represent the true 
schedule actually used by the engine/FCU. 
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FIGURE 5.3-3 

YAV-88 SHORT PERIOO MODE DAMPING RATIO 
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FIGURE 5.3-4 

YAV-8B PULL-UP IN WINGBORNE FLIGHT 

Gross Weight = 16.810 lb N F = 80% 0j = O # 6 F = 5° 
Mach 0.74 Altitude = 7,000 ft 
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Flight Data Altitude 

Symbol 

O 3,000 ft • 10,000 ft except as noted 

□ 25* flap 

Predicted 

7,500 ft, 1R.158 lb 
— — — Nonlinear moaei 
— — 1 Linear model 



Mach Number 
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FIGURE 5.3-5 

YAV-SB WINGBORNE 1 g ROLL PERFORMANCE 

3.2 in. Lateral Stick Input 
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FIGURE 5.3*6 
YAV-8B 

WINGBONE STEADY HEADING SIDESLIP CHARACTERISTICS 
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FIGURE 5.3*8 

Y/.V*8B RESPONSE TO RUDDER INPUT iN WINGBOP.NE FLIGHT 

Orjss Weight — 19.2o0 !b N F = 77% ^ = 3* 6 f = 5° 

Mucn 0 697 Altitude = 35.200 ft 
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FIGURE 5.4-1 
YAV-SB WET VTO 


5-28 















Time • sec 


QP21-0M244 

FIGURE 5.4-2 (Cont.) 

YAV-8B DRY VTO 

Nozzle Angle = 77* 
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The simulator results curing the dry VTO, Figure 5.4-2, appear to have 
a faster response than the actual engine. Ihis discrepancy is due to a 
calibration error in the throttle position of the flight aata. As a result, 
the initial engine fan speeds are offset. 

5.5 CONTROL SURFACE ACTUATOR/SERIES SERVO VA1IDA. ON 

A comparison of the frequency response characteristics of the control 
surface actuators and series servos is presented in Figures 5.5-1 to 5.5-4. 

The test data were obtained for small amplitude inputs while the predicted 
characteristics used the transfer functions presented in Section 2.5. Good 
agreement is shown for frequencies less than 1 Hertz, which are most pertinent 
to handling qualities analysis. Adequate gain margin and phase margins are 
shown. 



FIGURE 5.5*1 

YAV-flB STABILAiOR ACTUATOR - SAS MODE 
Small Amplitude Frequency Response 
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FIGURE 5.5*3 

YAV-8B AILERON ACTUATOR - SAS MOOE 
Small Amplitude Frequency Response 
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FIGURE 5.5*4 
YAV-8B YAW RCV SERVO 
S'Dall Amplitude Frequency Response 
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6 . PARAMETER ES T IMATION 

The importance of determining the aircraft stability ard control deri- 
vatives from flight test data has been recognized for many years. The flight 
determined derivatives are of much value to the flight test analyst. These 
results may be compared to other sources, such as wind tunnel test and 
theoretical calculations, to substantiate the predicted behavior of the air- 
craft. The extraction of stability ard. control derivatives from flight test 
data is a problem in parameter estimation. 

Over the past 30 years, considerable work has been done in extracting 
.he stability and control derivatives from flight test data. Various tech- 
niques have been investigated, some of which lend themselves to automation. 

As documented in References 2 and 3, techniques used by NASA include the 
Maximum Likelihood Parameter Estimation technique. To date, these efforts 
have concentrated on the conventional flight regime; however, with the 
acquisition c* the YAY-8B, NASA intends to extend this investigation into che 
V/STOL flight regime. 

This section contains suggestions on the forms of mathematical models to 
use in the estimation of V/STOL aircraft stability and control derivatives. 

In addition, some significant nonlinear effects, included and documented in 
the six degree of freedom mathematical model presented in Section 3, are 
discussed. 

6.1 JET BORNE AND SEMI-JET BORNE 

As discussed in Section 4, MCAIR uses the linearized representations of 
the aircaft in control system design and handling qualities analysis in the 
V/STOL flight regime. The comparisons to flight data of Section 5 show that 
the linearized moael is represe. ative of the actual aircraft. This form of 
the model was also used by CALSPAN in a study that assessed the feasibility 
of simulating the AV-8A in the terminal operation area, with the X-22A variable 
stability aircraft (Reference 4). A follow-on study that included flight test 
evaluation by 2 qualified AV-8A pilots concluded that the simulation was 
representative of the actual aircraft* 

In hover and at very low speeds (<35 RTAS) the simplifying assumptions 
detailed in the MIL-F-8330 Background Information and Users Guide (BIUG) , 
Reference 5, may be applied to the linearized equations of motion. These 
assumptions include: 

o hover (V T ® 0) 

o Vertical motion decoupled from pitching and horizontal translational 
motion 

o Yawing motion decoupled from rolling and lateral translational 
motions 

c Constant thrust magnitude and direction 
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Applying these assumptions to the linearized equations of Section 4, and 
recognizing that angie-of-attack, =, and sideslip angle, S, are no longer 
defined in hover, results in the simplified equations shown in Figure o.l-l. 

This form of the equations was used to analyze the flight test and flight 
simulation data available when formulating the hover and low speed dynamic 
response requirements of Reference 6. 

No matter which form of the equations is used, the following first order 
effects must be included in the formation of YAV-8B stability derivatives 
taken at fixed operating points in the V/STOL flight regime. Note that if 
the nonlinear model of Section 3 is used to establish a trim point and to 
formulate the stability derivatives, the following effects will be taken 
into account. 

o Thrust level and control effectiveness At low speeds the aero- 
dynamic controls are naturally ineffective and the primary source of 
control is the Reaction Control System (RCS). The effectiveness of 
thij system is proportional to engine RPM as shown in Figure 6 1-2. 

o Simultaneous demand and control effectiveness - Tne effect of multi- 
control usage is to reduce the effectiveness of the individual 
controls. This effect is shown in Figures 6.1-3 and 6.1-4 for the 
pitch RCS. 

o Thrust vector application point - The application point of the thrust 
vector varies with engine RPM as shown in Figure 6.1-5. 

One additional effect may be considered when evaluating the effectiveness 
of the RCS. AV-8B static tests using the individual Reaction Control ^alves 
(RCV ) 1 showed the direction of the RCV thrust vector to be a function of 
the valve opening. These AV-8B c a have been combined with TAV-8B valve 
orientation geometry to produce the estimated variations shown in Figure 
6.1-6. These variations are not included in the models defined in 
this report. Both the nonlinear and linear mathematical models u i the 
"design" RCV thrust vector angles. 

6 . 2 WINGBORNE FLIGHT 


The linearized equations of motion presented in Section 4 are applica- 
ble to wingbome flight. The effects of Vectoring ir Forward Flight (VTFF) 
should be invest-gated using stability derivatives based on perturbations from 
a trim point obtained using the nonlinear mathematical model. 
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FIGURE 6.1-1 

SIMPLIFIED LINEAR EQUATIONS OF MOTION 
HQVErx 
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FIGURE 6.1-3 

MAXIMUM RCS FRONT PITCH VALVc THRUST 

Normal Lift Dry Full Nose Up Control Demand 
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FIGURE 6.1-4 

MAXIMUM RCS TAIL PITCH VALVE THRUS^ 

N0RMU LIFT DRY FULl NOSE DOWN CONTROL DEMAND 
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7. CONCLUDING REMARKS 

A six degree of freedom FORTRAN batch simulation of a nonlinear mathe- 
matical model of the YAY-8B aircraft has been produced and documented. The 
documentation includes top level flow charts, a discussion of program struc- 
ture, subroutine interfaces, modeling equations, data format, a user's guide 
and plots of the over 17,000 data points used by the program. The model is 
based on the YAV-8B model used on the MCAIR manned flight simulator. The 
nonlinear model has been shown to be representative of the actual aircraft 
by comparison to both static and dynamic YAV-8B flight test data. 

Two simplified models have been suggested for use in parameter estima- 
tion, one of which is the traditional longitudinal and lateral-directional 3 
degree of freedom linearized model used in the analysis of conventional air- 
craft. A second model, applicable for hover and low speed analysis, is 
suggested. Aircraft stability derivatives, which characterize the aircraft 
throughout its flight envelope, are provded at 24 flight conditions. These 
derivatives were formulated by perturbing the nonlinear model after estab- 
lishing control settings corresponding to a 3 degree of freedom 1 'g T wings 
level trim. These stability derivatives include the effects of aerodynamics, 
inlet momentum and gross thrust. 

For reference a description of the major aircraft systems modeled in the 
FORTRAN program is included. 
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9. SYMBOLS AND ABBREVIATIONS 


SYMBOL 


DESCRIPTION 


UNITS 


ACU 

AND 

ANL 

ANR 

ANU 

AOA 

AR 

Autostab 

^Yacc 

b 

B. L. 

BOV 

C 

z 

Ca 

Cd 

C L 

C 1 

Cm 

C N 

C n 

C Y 

C. G. 

CPL 

Ctr 

CR 

CT 

DT 

e l to e 4 
EGT 

F A 

FCU 

FJI 

f N 

FOD 

F? 

FRL 

F.S. 

Fwd 

g 

GE 

GTS 

GW 

h 

HP 

HSA 

HT 

IGV 

Inbd 


Acceleration control unit 
Aircraft nose down 
Aircraft nose left 
Aircraft nose right 

Aircraft nose up - 

Angle of Attack deg 

Aspect ratio - 

Autostabilization or SAS 

Lateral acceleration at the accelerometer location ft/sec^ 


Wing span 
Butt Line 
Blow off valve 
Chord 

Mean aerodynamic chord 

Coefficient of axial force (positive forward) 
Coefficient of drag (positive aft) 

Coefficient of lift (positive up) 

Coefficient of rolling moment (positive RWD) 

Coefficient of pitching moment (positive ANU) 

Coefficient of normal force (positive up) 

Coefficient of yawing moment (positive ANR) 

Coefficient of side force (positive out the right wing) 

Center of gravity 

Centerline 

Compressor pressure limiter 

Center 

Root chord 

Tip chord 

Iteration time 

Quaternions 

Exhaust gas temperature. Same as JPT 
Force Axial (positive forward) 

Fuel control unit 
Flap jet impingement 
Force normal 
Foreign object damage 
Flat plate 

Fuselage reference line. Same as W.L. 

Fuselage Station 
Forward 

Acceleration due to gravity 

Ground effect 

Gas turbine starter 

Gross weight 

Altitude 

High pressure 

Hawker Siddeley Aviation 

Horizontal tail 

Inlet guide vane 

Inboard 


ft 

in 

in 

in 


Z Z , in 


in 

in 

sec 

°C, °R 
lb 


lb 


in 

32.174 ft/sec 2 


lb 

ft 
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SYMBOL 


Ijcc 

*xz 

1 yy 


k ay 

kcas 

KIAS 


Kp 

Kq 


Kr 

KIAS 

L 

LEMAC 

LIDS 

LP 

lswt 

Lt 

LVDT 

M 


MAC 

MCAIR 

MFC? 

mswi 

*v 

NED 

nf 

% 

NLD 

Ny 

0 

Outbd 

OWE 


?3 

*13 

PLA 

PM 

PEL 

P RCS 

PSA 


psia 

psig 


9. SYMBOLS AND ABBREVIATIONS (CONTINUED) 
DESCRIPTION 


Roll moment of inertia 

Product of inertia 

Pitch moment of inertia 

Yaw moment of inertia 

Jet pipe temperature. Same as EGT 

Lateral acceleration feedback gain 

Knots calibrated airspeed 

Knots indicated airspeed 

Roll rate feedback gain 

Pitch rate feedback gain 

Yaw rs~e feedback gain 

Knots true airspeed 

Left, Latitude 

Leading edge of the mean aerodynamic chord 
Lift improvement device system 
Low pressure 

Low speed wind tunnel (located at MCAIR) 

Left 

Linear variable differential transformer 
Mach number 
Mass flow 

Mean aerodynamic chord 
McDonnell Aircraft Company 

Mass flow calibration facility (located at M C A I R ) 
Minispeed wind tunnel (located at MCAIR) 

Rolling moment (positive is RWD) 

Pitching moment (positive is ANU) 

Yawing moment (positive is ANR) 

North, East, Down 
Fan speed 

High pressure compressor RPM 

Normal lift dry 

Lateral load factor 

Normal load factor 

Origin 

Outboard 

Operating weight empty 
Roll rate 

Compressor discharge pressure 
Fan delivery pressure 
Power lever angle 
Pitching moment (positive ANU) 

Pressure ratio limiter 
Reaction control pressure 
Power spindle angle 
Absolute pressure 
Gauge pressure 


UNITS 


slug - ft^ 
slug - ft2 
slug - ft^ 
slug - ft2 
•C, 

deg/ft/sec^ 

knots 

knots 

deg/ deg/ sec 
deg/ deg/ sec 
deg/ deg/ sec 
knots 
-» deg 
in 


lb®/ sec 
in 


ft - lb 
ft - lb 
ft - lb 

Z 

z 


lb 

deg/ sec 
lb/in2 
lb /in* 

ft - lb 

lb /in 2 

lb/in2 
lb /in 2 
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SYMBOL 


q 

Q. q 

q- 

R 

r 

RCS 

RCV 

r earth 

RM 

RPM 

RR 

RT 

RWD 

RWD 

S 

SAS 

SF 

SFO 

SLD 

SLW 

STO 

T,t 

*2 

T 3 

TED 

TEL 

TER 

TEV 

u 

V 

Veq 
VI FF 

V ? 

V/STOL 

V T 

VTO 

v 


w F 

WFT 

W.L. 

W N 

W.R.P. 

WOW 

WT 


9 . SYMBOLS AND ABBREVIATIONS ( CONTINUED) 
DESCRIPTION 

Pitch rate 
Dynamic pressure 
Freestream dynamic pressure 
Jet dynamic pressure 
Radius , Right 
Yaw rate 

Reaction control system 
Reaction control valve 
Radius of the Earth 

Rolling moment (positive RWD) 

Revolutions per minute 
Rolls-Royce Limited 
Right 

Right wing down 
Right wing up 

Wing area, Laplace transform 

Stability augmentation sv stem 

Side Force (positive out the right wing) 

Simulated flame out 

Short lift dry 

Short lift wet 

Short takeoff 

Temperature , thickness 

Temperature at engine face 

Compressor discharge temperature 

Trailing edge down 

Trailing edge left 

Trailing edge right 

Trailing edge up 

Velocity in X direction (positive forward) 

Velocity in Y direction (p ositive out the right wing) 
Equivalent velocity ratio 
Vectoring in forward flight 
Speed of sound 

Vertical/Short Takeoff and Landing 
Total velocity 
Vertical takeoff 

Velocity in 2 direction (positive down) 

Wind velocity, down comjr'ment 
Wind velocity, east component 
Fuel flow 
Fuel flow 
Waterline 

Wind velocity, north component 

Wing reference plane (parallel to W.L.) 

Weight on wheels 
Weight 


UNITS 

deg/sec 
lb /in 2 
lb/in^ 

lb /in 2 
in, - 
deg/sec 


ft 

ft- lb 


ft 2 , - 
lb 


°C, in 
°C, °K 
°C, °K 


f t/sec 
ft /sec 


ft/sec 

ft/ sec 

f t/sec 
ft/ sec 
f t/9ec 
lb/min 
lb/min 
in 

ft/ sec 


lb 
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9 . SYMBOLS AND ABB REYIAH. ONS ( CONTUSED) 

SYMBOL DESCRIPTION 

Xcg Fuselage station of c.g. 

YM Yawing moment (positive ANR) 

Z Z transform operator 


ONUS 

in 

ft-lb 


(SEEK SYMBOLS 


a 

» 

a 

6 

r> 

1 0 

Y 

A 

*A 

h 

6 h 

5 r 

°o 

5., 

id 

; «P 

n 

e 

e j 

^LE 

X 

P 

T 

♦ 

ii 

Q 

“d 

“sp 


Angle of attack 

Angle of attack rate 

Sideslip angle 

Initial flight path angle 

Flight path angle 

Incremental quantity 

Aileron deflection (positive is TED) 

Flap deflection (positive is TED) 

Horizontal tail (stabilator) deflection (positive 
is TED) 

Rudder deflection (positive TEL) 

Longitudinal stick deflection 

Lateral stick deflection 

Rudder pedal deflection 

Dutch roll damping ratio 

Short period damping ratio 

Fraction of semi-span 

Pitch angle 

Nozzle deflection 

Sweep angle of lea di ng edge 

Longitude 

Density of air 

Time constant 

Roll angle 

Yaw angle 

Rate of rotation 

Dutch roll natural frequency 

Short period natural frequency 


deg 

deg/ sec 
deg 
deg 
deg 

deg 

deg 

deg 

deg 

in 

in 

in 


deg 

deg 

deg 

deg 

slug/ f t J 
sec 
deg 
deg 

deg/ sec 
rad/ sec 
rad /sec 
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CONTRACT STATEMENT OF WORK 
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1. 0 INTRODUCTION 

As part of its V/STOL research program, NASA intends to conduct flight 
investigations of the stability, control and handling qualities of highly 
augmented V/STOL aircraft. Specific plans include the flight tests of a 
YAV-8B aircraft modified to include an advanced avijnics and flight control 
system for improved flying qualities and performance. 

As an initial phase to this program, NASA will conduct flight tests of 
the YAV-8B vehicle in order to extract aerodynamic and propulsion character- 
istics, update existing simulation models, validate handling qualities and 
design criteria, and to improve V/STOL flight test techniques. This program 
will also include tests using a static test stand location at the Drvden 
Flight Research Center, where flight tests of the YAV-8B will take place. 

In order to perform high quality parameter estimation and analysis of 
the YAV-8B characteristics, it is necessary to construct mathematical models 
of varying complexity and linearity from existing wind tunnel and flight test 
data. This procurement is intended to produce models and to compile existing 
data for the Harrier aircraft. 

2.0 SCOPE 


The contractor will assess existing wind tunnel, analytic, and flight 
data for V/STOL aircraft especially in VTOL-unique flight regimes in order to 
compile aerodynamic, propulsion, and system models. Detailed models of 
varying complexity will be derived for the YAV-8B aircraft for use in para- 
meter estimation as well as linear analysis programs. A simplified batch 
simulation of the YAV-8B, utilizing the best nonlinear representation of the 
airplane will be generated. Finally, flight data and wind tunnel predictions 
will be used to compare YAV-8B characteristics with flying qualities criteria 
and design quidelines. 

3.0 CONTRACTOR TASKS 


The contractor shall: 

3.1 Assemble and compile sufficient available Harrier aerodynamic, pro- 
pulsion system, and flight control system data, as well as available flight 
test data to support completion of tasks 3,2 through 3.8. 

3.2 Formulate the most complete model of the YAV-8B from the results 

of 3.1, to characterize the vehicle in hover, transition, cruise, and vectoring 
in forward flight (VIFF) flight conditions. 

3.3 Produce a simplified FORTRAN batch simulation of the YAV-8B 
(including engine) using non linear equations of motion, capable of execution 
on the Dryden Cyber 73-28 computer. 
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3.4 Formulate linear models of the YAV-8B for distinct flight conditions 
in hover, transition, cruise, and VTFF modes of flight. Select a sufficient 
number of flight conditions to characterize the aircraft throughout its flight 
envelope. 

^.5 Validate the models of 3.2 and 3.4 (with the exception of the VIFF 
mode), and the simulation program of 3.3 by generating time history responses 
and comparing them with actual flight data. Flake adjustments to the models to 
correct gross discrepancies. 

3.6 Recommend simplified (not necessarily linear) models appropriate 
for use in parameter estimation programs for each of the unique flight modes 
of 3.4, 

3.7 Provide MCAIR Report MDC A4637, Revision C, dated 11 January 1980, 
"YAV-8B Aerodynamic Stability and Control and Plying Qualities Report" at no 
cost. 


3.8 Provide MCAIR Report MDC A5032 Addendum No. 1 dated July 1981 
"YAV-8B Bulletin No. 158394" at no cost. 

4.0 REPORTING 


The contractor shall provide the following document^gpbn : 


4.1 Tables, functional plots, and equations of motion for the non linear 
model of 3.2. 


4.2 Linear coefficients and equations of motion for the models of 3.4. 

4.3 Listing and tape of the FORTRAN program of 3.3. Also, a complete 
description of the software, including top-level flow charts, program struc- 
ture, subroutine interfaces, and data format. 

4.4 Tine histories from non-linear and linear models for each of the 
flight r jnditions analyzed. 

4.5 A Final Report in contractors format which documents all analysis 
and results of this procurement. 

4.6 Monthly letter progress reports. 
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